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LAME f S WAVE FUNCTIONS OF THE ELLIPSOID OF REVOLUTION* 

By J. Meixner 


1. INTRODUCTION 


Lame's -wave f motions result "by separation of the wave equation 
in elliptic coordinates and by integration of the ordinary differen- 
tial equations thus originating. They are a generalization of Lame'B 
potential functions which originate in the same manner from the 
potential equation. Lame'B wave functions are applied for boundary 
value problems of the wave equation for regions of space bounded by 
surfaces of a system of confocal ellipsoids and hyperboloids. 

For general elliptic coordinates Lame's wave functions .have not 
been fully calculated so far. Except for a few general properties, 
not much is known about them. More consideration was given to Lame ' s 
wave functions for the case of rotationally symmetrical elliptic 
coordinates (called for short. Lame's wave functions of the ellipsoid 
of revolution) . However, even for these functions few resultB are in 
existence compared with those for the better known special functions 
of mathematical physics, such as cylindrical and spherical f unctions . 

The first more detailed investigation of Lame's wave functions 
of the ellipsoid of revolution was made by Niven (reference 1) who 
with their aid treated a heat- conduction problem in the ellipsoid 
of revolution. However, the numerical values of the coefficients 
of his series developments in terms of spherical and cylindrical 
functions as they are given, for the lowest indices contain several 
errors which were taken over into the report by Strutt (reference 2) . 
A more extensive investigation with a greater number of applications 
was xoade by Maclaurln (reference 3) • MSgllch (reference 4), whose 
mathematical investigation of Lame's wave equation is based on 
certain linear homogeneous integral equations, obtained results of a 


. Die Lameschen Wellenfunktionen des Drehellipsoids ." Zentrale 
fur wissenschaftliches Berichtsweeen der Luftf ahrtforschung des 
Generalluftzeugmeisters (ZWB) Berlin- Adlershof , Forechungsbericht 
Nr. 1952 1 June 1944. 
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more general character. Strutt (reference 2) gives a survey of the 
state of the theory of Lame’s wave functions in 1932# he also demon- 
strates on a large number of examples from acoustics, electrodynamics, 
optics, wave mechanics, and theory of wave filters, the manifold 
possibilities of application for these functions. 

Of the treatises published in the meantime, an investigation by 
Hanson (reference 5), which contains several new details, should be 
mentioned, aB well as a treatise by Morse (reference 6) on addition 
theorems, that is, on the development of tho plane wave and the 
spherical wave in terms of Lame's wave functions, furthermore, a 
number of treatises on the wave -mechanical treatment of the ion 
of the hydrogen molecule (reference 7) • Kotani (reference 8) 
deals with integral equations for .Lame’s wave functions. In 
particular, a treatise by 'Thu and Stratton (reference 9) should be 
pointed out which settles exhaustively the problem (treated so far 
only incompletely) of the continuation of the solutions of 
equation (2.4g) for large and small argument and shows in detail 
how the entire theory of Mathieu’3 functions results as a special 
and boundary case from the general theory of Lame 's wave functions . 
Finally, a treatise by Bouwkamp (reference 10) on the theoretical 
and numerical treatment of diffraction on a circular aperture is 
to be mentioned which, for the first time, contains more detailed 
numerical material concerning Lame's wave functions of the ellipsoid 
of revolution. 

The main task of the present report on Lame's wave functions 
of the ellipsoid of revolution will be to compile their most 
important properties in such a manner that these functions take on 
a form which facilitates their application. In this connection an 
investigation of the solutions of the ordinary homogeneous linear 
differential equations of the second order, which originate with 
separation of the wave equation in rotationally symmetrical elliptic 
coordinates, is of importance) further, it has to be determined what 
is to be understood in these solutions by functions of the first 
and second kind, their normalisation as well as the description of 
the behavior of these solutions in different domains of the 
independent variables, in particular, their asymptotic behavior. 

Here belongs also the indication of. a method of numerical calculation 
of these functions and the presentation of numerical tables. 

For the purpose of clarity it was necessary to generalize 
and supplement the existing material in some respects and to simplify 
some of the calculations and proofs. Therewith the theory of Lane's 
wave functions of the ellipsoid of revolution as a whole would seem 
to have reached a development equivalent to the theory of Mathieu's 


functions , and it probably even is somewhat simpler; Mathieu's functions, namely, represent not 
a regular hut a singular special case of Lame's wave functions. 


2 . TEE BASIC EQUATION 

2 .1 Rotationally Symmetrical Elliptic Coordinates 


In dealing with the rotationally symmetrical elliptic coordinates one must distinguish 
between those with oh late and those with prolate ellipsoids of revolution. Accordingly, one 
adds in the numbers denoting the formulas, the letter a or g, respectively, where both 
cases appear, x, y, z are the Cartesian, i , n, cp the elliptic coordinates, c is a 
positive constant. Then it follows that: 


x = c J(l - ti 2 )(I 2 + 1) cos 9 , y “ C \j (l - Tj 2 ) (| 2 + 1) Sin 9 , z 
0 = i < co, -1 i n = 1, 0=9^2* 



(2.1a) 


x = c \J( 1 - T] 2 ) (| 2 - 1) cos 9, y = C\l (1 - ti 2 )(| 2 - 1) sin 9, z = cItj 
1 ^ 6 < co^ -1 r \ < 1 , 0 ^ 9 = 2jt 


r ( 2 -ig) 


The symbols are the same as in the collection of formulas by Magnus and Qberhettinger (reference ll) 
to which reference is made also with respect to the transformation of the wave equation to 
elliptic coordinates and with respect to the separation of the wave equation. 


co 
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2 .2 Separation of the Wave Equation 

Solutions of the wave equation in three dimensions are to "be 
determined, (k = wave number.) 

Au + l^u = 0 (£.2) 

of the form 

u = f 1 (l)f 2 (q)f 3 (cp) (2.3) 

Then the ordinary differential equations 



(2.4a) 



are valid for f-^, ± 2 , and f q , . 

\ and y. are the separation parameters, first, they are 
assumed to be any complex numbers • They can only be determined 
for a given boundary value problem. In particular, y need not 


MCA TM Wo. 1224 


5 


"be an integer; this can "be recognized, for instance, in the treat- 
ment of an inside space problem in a sector 0 = <p Is cp Q of an 
ellipsoid of r evolution . 


2.3 Reduction to a Differential Equation 

The differential equation (2.4g) is designated as the basic 
equation. (2.5g) is identical with it; (2.4a) is transformed’ into 
it when I ’is replaced by fit and k^c 2 by -li^c 2 . Therewith 
the investigation of the differential equations (2.4a) , (2.5a), 
and (2.5g) is reduced to that of the differential equation (2.4g). 
The basic domain, however, is not .the same for all cases; it 
extends from -1 to 1 in the cases (2.5a) and (2.5g), from 1 to «> 
in the case (2 .4g) , whereas the basic domain of the differential 
equation (2.4a) in the transformation to (2.4g) will be changed 
to the domain from 0 to i<» (01- r^lse -ioo). It proves, therefore, 
to be necessary to investigate the differential equation (2.4g) 
in the entire complex I -plane. 


2.4 Transformations of the Basic Equation 

The basic equation represents a special case of the linear 
homogeneous differential equation of the second order with four 
extra essential singularities, two of which are made to join to 
one essential singularity. The latter is at infinity, the two 
remaining extra essential singularities are at 1 and -1. The 
present investigation of the basic equation will start with 
connecting its solutions with o the solutions of limiting cases of 
the basic equation. For Ire'' = 0, the basic equation is trans- 
formed into the differential equation of the spherical functions 
and their associated functions or, as they will be called here, 
of the general spherical functions. If one lets the two singu- 
larities at 1 and -1 combine into a single singularity at 1=0, 
there originates, aside from an elementary transformation, the 
differential equation of the cylindrical functions. This is 
brought about by the substitution 


i 


(2-7) 


and the abbreviation 


7 = ltc 


( 2 . 8 ) 
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if one then performs the limiting process y ->0. From (2.4g) there 
originates with ■ 


fl = (t 2 * - 7 2 y/ 2 ^ vft) (2.9) 


the differential equation 

v = 0 (2.10) 

In the transition from (2 .4g) to (2 ,4a) £ is transformed, into 

itself and y'~ need only he replaced by -y^ . For large distances, 

p o p p ~p 

that is, r = sr + y + z » 



^2 = k 2 r 2 


/'■vr2 + v 2 \ r 2 

1 - + 


(2.11a) 


£ 2 = -k 2 r 2 


i + c/sji\ 

r 4 Vr 4 / 


(2.11g) 


are valid. 


Another important- limiting case of the basic equation occurs 
if, of the two singularities of the basic equation located at finite 
distance, one or both move to infinity. Then the differential 
equation of Laguerre's and Hermite's orthogonal functions, respec- 
tively, is formed. This limiting case will yield, the asymptotics 
of the eigenvalues and eigenfunctions for large absolute value of y . 


2 ,s Connection with Mathieu's Functions 

Mathieu's functions are, in connection with Lame's wave functions, 

obtained in two ways. They appear, as is well known, in the separa- 
tion of the wave equation In the coordinates of the elliptic cylinder 
and must, therefore, also appear in the limiting case of Lame's wave 
functions for the ellipsoid with three axes when one axis becomes 
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infinitely long. However, Mathieu ' s differential equation is also 
obtained, except for an elementary transformation, if p in (2.4g) 
is set equal to tl/2 . This also indicates that it is useful, to 
consider the basic equation not only for p that are integers, 
but rather for arbitrary coefficients V and p . The theory of 
Mathieu's functions is, therefore, a special case of the theory of 
Lame's wave functions of the ellipsoid of revolution. Although 
the present report does not yield new results of Mathieu's functions, 
it demonstrates how they fit into a more general picture.. 


3 . SPHERICAL AWD CYLINDRICAL FUNCTIONS 
3*1 A Few Formulas for Spherical Functions 


The most important formulas and theorems for spherical and 
cylindric fimctions needed below are compiled and a few estimates 
for these functions are given, which will be necessary for con- 
siderations on uniform convergence of certain series in terms of 
such functions. Magnus and Oberhettinger (reference H) ia again 
referred to concerning the notation and additional formulas. The 

general spherical functions P^(|) and Q^(!) both satisfy the 

differential equation 


.4. 

di- 


(1 




v(v + 1) 



p£(!) 


= o 


(3-1) 


and both satisfy the recursion formula 


(2V + l)gpj(0 = (V - p + 1)P V+1 (!) + (v + n)Ejf -;L (6) (3.2) 


from which by three times repeated application 


I 2 ^(I) 


& 


2 ) (v - p + l) p. , . 

(2V + 1)(2V + 3 ) — V+2 ' 


2 V‘ 


2V - 2u c 


(2V - l)(2v + 3) 


^(D 


Av + u)(v + u - l) r u ft\ 

(2V - l)(2v + 1 ) —V -2 ' 


(3-3) 



is obtained. Further, the series presentation which converges for ]| j > 1 is given: 


co 


£< S > = * 1/2 e 1 - 1 * 1 ( t 2 - lf /2 r * ( 2|)' V ’ 1 \ + ?S . + . 1 ) ( g |) - 2 p 

Z— . r(v + p + |)r(p + i) 


(3 A) 


p 

where arg (I - 1) = 0 when I is real and >1, arg 1=0 when I is real and >0. In order 
to obtain the uniqueness of the general spherical functions a branch cut is put from £ = - w 
over -1 to 1. 


3 *2 A Few Estimates for Spherical Functions 
For the general spherical functions of the second kind the integral presentation 


2^r (u + iVZ /id+,0+,1-,0-) 

i (n-v)ni v . 2/ \G/ /g "'j p/2 -v-ii -1 

V 5) “ 16 4* sin (V + V Z 


J 


U V+ “ (1 - 


X (l - i) 
V Z / 


-1/2-H 


du 


(3-5) 


is valid. 
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f—g Vs, 

Therein z = ? + y| - 1; the sign of the root must he selected so that jzj = 1 
(one need hardly he afraid that z could he mistaken for the Cartesian coordinate z). la: 

order to estimate (^ (£) for r-^too one forms the absolute values of the individual 

factors ahead of and in the integral. !uj r in the integral is replaced by its maximum value. 
This latter can, by suitable selection of the path of integration, be made smalle r than (l + S) r 

for r > 0, smaller than S r for r < 0, where 8 is an arbitrarily small but fixed positive 
number independent of r, £ . Then one obtains 



r (i + s) r i 

-j r ‘ < ( s> 

L 8 J 


for r = 0, 1, 2, . . . 
for r = -1, -2, . . . 


(3-6) 


$^(|) is a restricted positive function of I independent of r in each closed domain 

excluding the points i = il,oo. Only the case of v + p being an integer requires special 
consideration. For V + p = 0, 1, 2, . . , it can be demonstrated that the estimate (3.6) is 

retained; for v + p = - 1 , -2, . . . the estimate above is valid fir (^)j r(v + r + p + l) 

instead of (£ ) • The case p = -i, -2., ... is not to be excluded, as can be shown easily. 


The general spherical functions of the first kind are given by the integral representation 




f)( l+,z 2 -) 

z - v -p-! I u V + p 

d 


(x - (1 - j )- 1 / 2 -* 1 au 


(3.7) 


VO 
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It is valid for -2* < arg z < 2jt with the provision that the path 
of integration for Jt < arg z ? < 2jt and -2n < arg z < -it, respec- 
tively , leads past the left of the point u = 0 and is also to he 
returned there . The estimate has to he made as above . The maximum 
of |uj r . is, if the path of integration is suitably selected, 
smaller than (Izj^ + &) r for positive r's and smaller than 5 r 
for negative r's, where 5 is a number of the conditions indicated 
above, so that 

(6) S (jzj + for r = 0, 1, 2, . . . 

v+r \ ! z I / v 

( 3 - 8 ) 

= jzj' r f^U) for' r = -1, -2, -3, . . . 

respectively, where \|^(£) and ^(?) are positive, restricted 
functions independent of r in each closed domain excluding the 
points I = ±l,oo . 

3.3 A Few Formulas for Cylindrical Functions 


For the following it is more convenient to introduce not the 
cylindrical functions themselves but rather the functions 


“ V 2 ? F v+l/2^^ ~ 'J 2 1, \ + l/2& (3*9) 


They both satisfy the differential equation 




1 


5 d£ 


v(v + 1) 



= 0 


(3.10) 
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and the recursion formulas 


ZL ± i = t 


2V + 1 




£ d£ “ 2V - 1 V " 2 (2 V - 1) (2v + 3 ) v 


. y . Tjf /o 
2V + 3 V+2 


Sv 


-i— i- ijr - \jr + l rS?V. JLJkl i|r + i \|f (-s ton 

2 *V 2v - 1 V-2 (2v - 1)(2V + 3 ) V 2V + 3 V+2 


Besides, the simpler recursion formula 


Pv 

j - Vi ( ^ + z v-i ( ^ 


(3.13) 


is to he noted for the cylindrical functions from which (3*12) is 
obtained by repeated application. Finally, the consistently convergent 
series development 



(3-14) 


is given with arg £ « 0 if $ is real and >0 .• 

3*4 An Estimate for Cylindrical Functions 

For the cylindrical functions Z (£) one obtains by repeated 

v-J-r : 

application of C 3 *13) ; 





\ r r (v + r) 

r - 1 ft) 2 , 

(r - 2)( r - 3) (tf 

) r (v) 

(v + r - l).v.i!\2/ (• 

v + r - l)(v + r - 2)v (v + 1)2 !\ 2 / 


z v (£) 


H 

ro 


AY' 1 r.(Y_t., r) 

W r (v + i) 


1 - 


r - 2 ft 

(v + r. - 1) (v + 1) 1 1 V 2 , 


(r - Djr - 4) 


(V + r - 1){V + r - 2)(v + l)(v + 2)2! V 2 


Vl^ 


(3.15) 


The number of the sum terms in brackets to be included is ~ + 1 for even and for 

d. c. 

odd positive r's. If Re v » 1, one obtains from (3*15) the following estimate 


W 0 


2 j r f(v r ) 

l\ r ( v) 


Zjt) F n {i!C!) + z v . 1 (0|F 1 (iiCI) 


(3 .16) 


Therein r = 1, 2, 3, • • • • A corresponding estimate may be obtained for r = -1, -2, - 3 , 
Since the singular points of the cylindrical functions lie at 0 and 00 , the function • 


z (t) 

v+r b 


r( v+ r) 


H 

ro 

ro 


(3.17) 
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is in every closed domain, excluding the points 0 and oo, 
restricted (considered as a function of t)} the upper limit does 
not depend on r» If Re V < 1, this .relation is valid at least 


for such r s for which 
is valid for negative 


r's, 


Re (v + r) = 2.. A corresponding estimate 



4 . THE X-IWCTIOlfS OP THE FIRST AND SECOND KIND 
4.1 Definition of the X-Functions of the First and Second Kind 


Since for y = 0 the "basic equation (2 ,4g) is joined to the 
differential equation (3»l) of the general spherical functions, it 
suggests itself to develop the solutions of the "basic equation in 
terms of spherical functions. One formulates the two at first 
formal series 


iV (y)P^ U) (4.i) 

v / v, r v+r 

r--co 


X^ (2 ) (|;r) « i r a^ (y)Q? (fi) (4.2) 

v b v,r v-i-r 


and attempts to determine the coefficients a 1-1 (7) and the index v 

v,r 

in such a manner that these two series formally satisfy the "basic 
equation and converge. The further problem will be to investigate 
the convergence properties of the two series (4.l) and (4.2) in 
order to determine that, for the two series, one has to deal with 
analytic f -unctions which, in general, are linearly independent 
solutions of the basic equation. 

For the coefficients (7) the indices v ju and the argument 7 

V ,v 

will be omitted where there is no danger of confusion; the same 
applies to the coefficients to be introduced later for series 
developments of a similar kind.. The summation index r assumes only 
even values. The term with r = 0 in the two series (4.1) and (4.2) 
is designated as the principal term of the series. In the solutions (4.1) 
and (4.2) of the basic equation an, arbitrary constant factor remains 
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f undetermined. It may "be determined In some way. Then the series (4.1) 
■is denoted as X-f unction of the first kind and the series (4.2) as 
X- function of the second kind with the argument 1 and with the 
indices v,|j. with the parameter 7. It will he found that the index v 
is determined "by the separation parameter X; more accurately, 


there exists a functional relation between 
which is expressed by 


X = X^(y) 
V 


X, 


(i, and 7 


(4.3) 


The series (4.1) and (4.2) are now inserted in the basic 
equation, the differential quotients of the spherical functions are 
eliminated by means of the differential equation of the spherical 
functions (3*1), and the factor of the spherical functions is 

eliminated by application of (3*3) • Then there appears an infinite 
sum of spherical functions with coefficients independent of | which 
is equal to zero. The disappearance of the individual coefficient is 
sufficient to this end. This leads to the conditional equations 


2®i^r ~ c -ir a r~2 + Pr a r+2 
7 


(r = 0, 1'2, 1"4, 


.) (4.4) 


with the abbreviation 


3> = -x + (V + r + 1)(V + r ) + 7 2 r ) z£ t£ ~ 1 

r (2V + 2r + 3)(2V + 2r - 1) 


(y+r+u + 2 )(y + r + u + l) 
P r “ (2 v + 2 r 5 ) ( 2 v + 2 r + 3 ) 


r 0^5) 


q = (v ± r “ u) (v + r - u - 1 ) 
r (2V + 2r - l) (2v 2r - ,3) 
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4.2 General Qualities of the Coefficients a r 

The recursion formula for the coefficients a r .is interpreted 

as a difference equation. In order to avoid complications, the 'case 
of real fractional values of one-half for V is completely excluded . 
and the case of real integers for v + M- and v - n, respectively, 
is postponed. Concerning the "behavior of the coefficients at 

infinity, a simple formulation can he obtained according to 
Kreuser (reference 12) . The equations 


lim sup v r j la - r 1 _ 2 
r— > oo v r i j 7 j 


(4.6) 


or 


lim sup 
r— 




(4.7) 


are valid. 

If the behavior of the coefficients a r at infinity is given 
by (4.6) at least for negative or- positive r's, they increase too 
strongly to make a convergence of series (4.1) and (4.2) possible. 
Therefore, a solution of the difference equation (4.4) is to be found 
which shows the behavior (4.7) for r— > oo as well as for r — *-«=. 
Althoxigh there always exists an exact solution which behaves 
for r--~? - oo as indicated in (4.7), this solution will in general 
exhibit for r — > - °° the behavior (4.6). Only for certain distinct 
values of the parameter V (free so far), the behavior (4-7) 
prevails for both r — 4 “ and r — inversely, in this manner 
distinct values of X are coordinated to each value of V. For 
7 = 0 the conditions are part icularly .simple . There becomes for 
all r's 



( v + r + 


l) (v + r) 


a r = 0 



•) 


Thus X can, for a given V, assume any of the values (V + r + 1) x 
(V + r). It is determined by the requirement that the series (4.1) 
and (4.2) should be reduced to the principal term for this case, 



1 6 
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which leads to X = v (v + l) . Now it is further required that 
under ^(7) always the value should he understood which goes over 

to v (v 1) for 7~3> 0. The existence of such a distinct 1- value 
to each given V , p, and e and- its uniqueness will not he proved, 
here j it follows from the method of calculation given in section 6 
for the determination of X . 

From here on. the coefficients a 1-1 (7) will always represent 

V,r 

that solution of the difference equation (4.4) which shows the 
behavior (4.7) for r — >i co, belongs to the value 1^(7), and 

therefore has the boundary values 


lim a r = 0 (r = 'i'2, 1'4, 1'6, . < .) (4.8) 

7=0 


Furthermore, the constant factor which is arbitrary in the coeffi 
cients a^, may be determined in a given manner. 


4.3 Convergence of the Series Developments of the X -Functions 
of the First and Second Kind. 

From the estimates (3.6) and. (1.8) as well es from the boundary 
values (4.7) for r —>t°° there follows immediately that the 
series (4.1) and (4.2) in each closed domain, which does not include 
the points £ = i'l, 00, will converge absolutely and uniformly. One 
may further conclude that the series (4.1) and. (4.2) will converge 
as well as the exponential scries. Since the individual terms of 
these scries are analytic functions in this domain, there follows 
from the uniform convergence that the sums of the series themselves 
will again he analytic functions, the singularities of which can. 
lie only at | = ±1,«, furthermore, that the series can be dif- 
ferentiated termwise, and therewith the fact that the functions 
represented, by these series are real solutions of the basic equation. 
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4.4 Further Solutions of the Basic Equation and Their Relation 
to the X-Functions of the First and Second Kind 

Between the general snherical functions P* 1 , P - ^ , p^ 1 . 

~v ~v --v-1 

P ^ , Q , Q 11 all of which satisfy the same 

~-V-l v V ~-V-l “-V-1 

differential equation there exist, in general, six linear relations 
independent of each other. They can be generalized for the X-functionB 
of the first and second kind. To this end several relations for the 

coefficients a^ (7) will be derived. 

v,r 

The system of equations (4.4) and the system of equations 
originating from it by the substitution V — > -v-1 and r—>-r are 
identical because of 


-v-l,-r 


_ .11 
= $ 


V,r 


-v-i.-r 


l-i 

- V 


v,r 


Due to the uniqueness of the solution there follows from it 


IT (7) — V (7) 

-V-i V 


(4-9) 


Furthermore, the constant factor which is arbitrary in the a r 's can 
be determined in such a manner that 


-v-l,-r 


(7) 




V,r 


(7) 


(4.10) 


The system of equations (4.4) and the system of equations 
originating from it by the substitution ijl— become identical 

if one introduces in the latter instead of the a 11 's the values 

v,r 


( 7 ) = L(v ±J L ± - tL ± J± r,-(v_iji_+.li . 

v,r r(v + r - u + 1) r(v + \x + 1) v,r ' 


(4.11) 
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They are determined so that b Q = a Q . Then the equations 


^ (r) » a ^ (y) 


(4.12) 


v,r 


v.r 


and 


x->) - ^y(r) 


(It. 13) 


are valid. 

After these preparations, at first a relation between X (| }y ) , 
X^^(|j 7 ), and ^ ( 7 ) is -derived, since according to (4.1) 

and (4.12) 


\ Al) k ») 


\ r jjl . -(a , . 

> i b (y)P (?) 
X v,r v+r 5 


(4.14) 


r=-» 


is valid. 

If one expresses in this equation the spherical function P ^ (i ) 

" V+r 

by F (| ) and Q (| ) (under consideration of (4.1), (4.2) 
v+r ~v+r 

and (4.11)) the required relation 


j _ r(v - u + 1) 
X V *’ 7) r(v + u + 1) 


,-M- (l) / \ 

(S ir) 


p -M-Tti . ' |J- ( 2 ) , 

sm pit (5 }y) 


ft 


(4.15) 


will .be found at once . 

In exactly the same way there result the formulas 


x^ (1) U;y) = x u(1 ) (e;r) 

-v-i V 


(4.16) 


Z? (S) (S;7) = e- 2 "” 1 C M-tL ± 1 1 y?fhi;r) 

r(v + U + 1 ) V 


(4.17) 


^(|jr) sin (v + n)s - sin (v - p)jt - Tte^ 1 cos Vit /^(g; r ) (4.3JB) 


J*(2), , . -'Vjti u (2) % 

(“1^7) - “ e (ijy) 


(4.19) 


d v (1) 


(-|;r> - ' ” :l ^ v'" 1 ( ' i '/ ) - | Bln (v + nine' 11 ’ 11 (tir) (4.20) 


For the last two relations the upper or lower sign is valid depending on whether the 
imaginary part of g is positive or negative . 


The relations (4.15) to (4.20) are identical with corresponding relations for the general 
spherical functions. This was obtained by the determination of the arbitrary constant factor 
in the a r 's. That determination is not yet unique j it leaves a certain latitude. It is 
compatible with (4.10), (4.11), and (4.12) if one sets the coefficient of the principal term 
equal to one ) but then it may happen for certain combinations of values v, p, and y that 
all a^'s with the exception of a Q become infinitely large. Other determinations still 


possible within this latitude would be the requirements that 


or 


\ 2v + 1 

2v + 2r 1 ar r = 


00 

\ 

\ 

/ 

/ 

p— - 00 



00 


\ 


£ 

30 
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4 .5 General Relations between the X-Functions 

The following general relations are valid for the general spherical functions in the 
case j | j > 1 


-V -V 


(4 . 21 ) 


P^(|e llti ) = e W3ti p£(|) 


2i 

Tt 


sin nl (v + i) 

\ 2 / 

t -| \ 

sin si lv + — } 
\ d / 


sin (v + jj.)jte" Cl a+ l/2^)«i 


( 4 . 22 ) 


l is therein an integer, either positive or negative. The coefficients in (4.21) and (4.22) 
remain unchanged if v is replaced by v + r, r being any even number. If one replaces v 
in (4.21) and (4.22) by V + r, multiplies by i r a r and forms the sum over all even r 
from - «> to co , there is formed because of (4.1) and (4.2) 


= e ”* 1 


jtl(l) ( 


i;y) - 


f 

0 . sin it 2 ( V + -i 

2l v cj_ 

71 sin 3Ti'v + 

■v 2/ 


sin 


(v+ |i)jte"" X^(^) (^ *y) (4.23) 


X^( 2 )( ? e^ i ; 7 ) = e -Z(v+l)«i X* a ( 2 ) (| ;-v) 

V V 


(4.24) 


MCA TM No. 1224 



MCA TM M. 1224 


21 


These two general relations which are valid for j g [ > 1, are 
in the case 2 = 1 transformed not to exactly (4.19) and (4 .20); 
Se* 1 -, namely, represents in ( 4 . 23 ) and (4.24) an increase of the 
argument hy it, whereby under certain conditions the branch cut 
may be passed, whereas the argument of -£ in (4.19), and (4.20) is 
obtained by choosing such a path from £ to - i that the branch 
cut extending from -<» to 1 will not be passed. 


5 . THE -Z-FUNCTI0NS OF THE FIRST TO FOURTH KIND 
5*1 Definition of the -Z -Functions of the First to Fourth Kind 


If the two extra -essential singularities of the basic equa- 
tion (2.4g) are made to join, as indicated in section 2, there 
originates, aside from an elementary transformation, Bessel's 
differential equation. It therefore suggests itself to attempt a 
solution of the basic equation also by series developments in 
terms of cylindrical functions . The functions defined by the series 
(which are at first formal) 



( 5 ^ 1 ) 


(5-2) 


are defined as E-functions of the first and second kind. In these 


series arg £ < it; 


arr 


Z-Z 2 . «- 

2 


if 


i 


2 p 

arg £ = arg y . 


Substitution of these series into, the basic equation (2.4g) (it is 
best to insert it into the transformed form (2 ,10) of the basic 
equation) , elimination of the first and second derivatives of the 
indices + v+r and n v+r by means of (3*10) and (3*11), and 

removal, of the denominator by means of (3*12) leads finally, 

exactly as in the X-functions of the first and second, kind, to a 

three-term recursion system for the b^ (y) . It agrees with the 

v ,r 

recursion system (4.4) for the a 11 (y), if the 

V,r 


H there is 



22 


MCA TM No. 1224 


replaced by -p. The solution of the recursion system differs from 

the indices b^ ( 7 ) defined in (4.11) only by a constant factor; 
v,r 

this factor is selected to equal one. Therefore the relation (4.3) 
formerly found between the index v of the generating functions 
in ( 5 . 1 ) and ( 5 * 2 ) and the separation parameter \ has to be 
assumed also in this case . 


As Z --functions of the third and fourth kind one defines 

< (3) (?;r) =< (1) (t;r) + < <2, (?;r) (5-3) 

af>( S ;y) = (i ;y) - tZy' 2 ' 1 (f jr) (5-*0 


They have the same relation to Henkel's functions as the 
Z-functions of the first and second kind to Bessel's and Neumann's 
functions . 


5-2 Convergence of the Series Developments of the Z-Functlons 
of the First and Second. Kind 

It must now be demonstrated that the series (5»l) and (5*2) 
converge -uniformly in a certain domain. One starts from the 
estimate ( 3 * 17 ) and from the boundary values (4.7) which are also 

valid for the b 11 . There results 

v,r 


lim sup 
r_> oo 


t/i* 


a tr z vdr 


< 




(5-5) 


The convergence is uniform. Thus the series (5*1) and (5*2) con- 
verge uniformly and absolutely in the entire domain ||| >1 with 
the exclusion of the infinitely distant point; they represent 
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therefore analytic functions, can he differentiated any number of 
times termwise with respect to and satisfy the basic equation. 
Only in special cases these series converge also for ||j =1. 


5.3 General Relations between the E-Functions 

The transition to. various function branches over the branch 
cut from through to 1 is made possible by the general relations. 

They can be obtained corresponding to the case of the X-functions 
from the general relations valid for the separate series terms, 
thus for the cylindrical functions. (Compare Magnus and Oberhettinger, 
elsewhere.) They read for )£| > jyj that is, j Ij > 1 


2 ^( 1 ) (e ln H;y) «' e Wlti ( 5 - 6 ) 


gg ( 2 ) (e^^i^ jy) _ e - Z(v+l)rti 


2 ie -lni /2 gin + cot (£ }y) ( 5 * 7 ) 


gM’CS) ( 0 Zrti£ *y) - 


_ e -Irti/2 


i /2 3ln (l - X) ( V + IV 


sin (v + gji t 


t (3, «;r) 


. sin Z(v + ijjt 

5 -Zjfi /2 e -(v+l/2)jti i 2 Z_2n(4)^. 7 j 


sin 


(5*8) 


z !j(' 0 ( o T.iti|- , 7 ) _ 


. a -*" 1 / 2 .( VI ft * 


sin (v + §) 


Jt 


, sin ( 2 + 1 ) (v + ~) « 

+ g-Ziri/2 _ ' 2/ 


sin |v + it 


zt (l ° ;r) 


(5.9) 



5. 4 Asymptotic Developments of the Z -Functions 

ro 

In order to obtain asymptotic developments of the Z-f unctions for l£j » l, it suggests ^ 

itself to insert the asymptotic developments of the functions t v+Y ,(0 ’and rjy + (£) into (5.1) 

and (p* 2 ) and then to interchange the summations. Since 


Vr ( ^ t 


v+r 



( 5 . 10 ) 


(-n < arg £ < 2jt and -2it < arg ? < it for the upper and lower sign, respectively) where , for 
abbreviation 


(v + r + — >v) es 

V 2 y plr(v + r - p + 1) 


( 5 - 11 ) 


there results in this manner 


< (3,1 ° (Or) - (f, 2 - r -)' 1 ' 2 ) _oM_ + o( El -«) 

P=5 ft2«) p 


'M-l 


(5-12> 


(-it < arg £ < 2jt and the upper sign for Z^3)(^ ‘ } y) } -2it < arg £ < it and the lower sign 

JJ ^ 1}. ^ ^ 

*7 ) ) • Ihe coefficients C(p) are defined by the absolutely convergent' series 


C(p) = ^ (v 


3fr=- 00 


r + 



i r b. 


M- 


v. r 


( 7 ) 


(5-13) 
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This derivation is not accurate since the asymptotic develop- 
ments (5 .10) are further dependent upon the condition )•£ j » |v + r | , 
and this condition is not satisfied for all series terms of equa- 
tions (5.1) and (5.2), respectively, since the sum has to he formed 
over al l r’s from -00 to » . The fact that the developments (5*12) 
are valid nevertheless is due to the behavior at infinity of the b r ’s 
(compare equation (4.7)) according to which the series terms with 
sufficiently large values of r do not contribute noticeably to * 
the 2-functions . 

Equation (5. 12) is proved as follows. 

According to general theorems on the asymptotic behavior of 
the solutions of homogeneous linear differential equations, the 
coefficients of which are polynomials, one obtains asymptotic series 
for the solutions by going into the differential equation (2 ,4g) 
with a formulation of the form (5*12) and attempts to satisfy it 
formally. This yields for the present case for the coefficients C(p) 
the four-term recursion system 


(p + l)C(p + 1) + 


(p + l)p + r-A 


c(p) + 4y a (n + p)c(p - 1 ) 


+ 47 2 (n + p)(n + p - i)c(p - 2 ) = 0 

C(-l) = C(~2) =0, p » 0, 1, 2, 3 , 


(5-14) 


from which they can be calculated recursively . This recursion 
system, however, is satisfied just then when the series (5*13) 
are substituted for the coefficients C(p) . This substitution 
leads after slight transformation to 


X i**. Jr) J r(v 

r=-oo 


V 7 * 


r(v+r-p+l) 


(v + r + 1) (v + r) + 7 ~ X. 


ro 

os 


r( v + i- +.p - i) i .. 2. 


r(v+r-p+3) 


+ p) 


( V + r + 1) ( v + r) + (p + 1) (p - 1) 


= 0 ( 5 - 15 ) 


for p » 0, 1, 2, 3, ... and these relations between the h r 's can he simply derived 
hy multiplying equation (4.4) hy i r r(v + r + p + l)/r(v + r - p + l) , hy forming the sum 
over r from -« to «, and finally replacing p hy -p and therewith a ^ r (r) *7 b v,r^ 7 ^ * 


Thus the asymptotic series (5.12) with the significance of the coefficients C(p) given 
in equation (5-13)" are actually asymptotic solutions of the basic equation (2. kg); it is now 

easily understood that they represent asymptotically those solutions ’ (£ir) • 


In a special case the developments (5-12) are even convergent: then namely, when v is a 

real integer = n. Then the series (5-10) are broken off j the functions ^^(0 and n n+r (0 

are elementary functions. Since the series (5-1) and (5*2) converge uniformly, the summations 
can he exchanged after substituting equation (5*10) in equations (5*1) and (5 * 2 ) j therefore 
the series (5*12) are convergent. Their domain of convergence is j | j > 1, that is, |£ { > jlj* 
One can also easily understand with the aid of equations (6.6) and (6.7) that, for p > n, 
the coefficients C(p) are of the order of magnitude y a , with a being the smaller one 
of the two even numbers among the four positive numbers p - n, p-n+ 1, p + n+ 1, 
and p + n + 2. For the case when v and p are positive integers and V > p > 0, the 
series (5.1) of the Z-f unctions of the first kind converge for all finite £ . 
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5*5 Further Solutions of the Basic Equation and Their Relation to 
the ^-Functions of the First and Second Kind 

It is immediately clear that with the Z-f unctions Z y ^ ( £;y) ' 

and Z the functions ( £;y) and Z~ v ^^ 1,2 \t,;y) also 

are solutions of the "basic equation. Since there exist only two 
linearily independent solutions of the Basic equation, it must "be 
possible to express all solutions linearily "by two of them. Because 
of the two relations 


^.. v = "cos v*n^(£) - sin vni|r (£) 
n -v-i^) - COS Vm|f v (£) - sinVjtn (£) 


f- (5*l6) 


there follows with the aid of equations 
definitions (5.1) and (5-2) 


= - sin vnZ 1 ^ 1 ^ (£;y) 
-v-1 V 

Z?^Uir7) = COS vn-z^^jr) - 

-v-i . v 


(4.12) and (4.10) from the 

- cos VitZ 1 "^ 2 ^;/) (5.17) 

v 

sin VrcZ^ 2 ^(£;r) (5 .18) 


In order to express the functions £^-(1,2 )(^. y ) £^(l)(^. y ) 

it will he practical to use the asymptotic series j 

it is sufficient to limit oneself to the first term of the series. 
Then there Becomes 


z ^(3,4) 


(Ur) 



iV 


v r 


(7) 


1 + o(|s| - 1 )] 


(5-19) 


l 
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The only difference for 

is that here a* 1 (7) takes 

v,r 

immediately that 


the asymptotic series for Z^^ ^ (£jx) 

the place of b^ (7) • There follows 
v,r 


z' u(1) (c;r) 21 

Y»=-rc 




(i .1, 2) (5.20) 


By combination of equations (5. 17), (5*18), and (5*20) finally, 
also the solutions z"^^(^7) of the basic equation can be reduced 

to the two solutions (t, }y) • As special cases of equations (5 

and (5*18) 


z^<Sir> - i®'”' 1 z^ (3, (Ur) 

-V-l V 


(5.21) 


should be noted. 


5 .6 Laurent -Developments for X- and E-Functions 

The X- and 2-functions were introduced wholly independent of 
each other. Since they all are, however, solutions of the same 
differential equation, it must be possible to express, for instance, 
the Z-functionB of the first to fourth kind in general linearily by 
the X-f unctions of the first and second kind. It will appear that 
in general the Z -functions of the first and second kind are not 
proportional to the X-functions of the first and the second kind, 
respectively. Thus, it is not possible to define simply functions 
of the first and second kind for the solutions of the basic equation; 
it must always be added whether one is dealing with X- or Z-functions 



The problem to express the Z-functions by the X-functions can be solved -partly by 
comparison of the general relations already known for both hinds of functions” (4.24), ( 4 . 25 ), 
(5-S), (5*7), (5*8), and (5-9) • A complete solution of this problem, however, is obtained 
in the following manner. The series (4.1) and (4.2) for the X-functions and equations (5.1). 
and (5.2) for the Z-functions have the common domain of convergence 1 < j|| < oa- They even 
converge uniformly in the domain bounded by the two circles |£| =1+5 and |i| -l + v, 
where 5 is an arbitrarily small and -3 an arbitrarily large positive number. Since the 
series terms of these developments are analytic functions regular in this domain, they can 
be developed in Laurent- series which converge in this entire domain. According to 
Weierstrass 's double series theorem one may exchange the summation over r and the 
summation of these Laurent- series and thus obtain the X-functions and the Z-f unctions in a 
representation by Laurent-series- From the comparison of these Laurent- series then results 
the representation of the Z-functions as linear combinations of the X-functions of the 
first and second kind. 

The performance of the transformations described Just now yields, using equations (3.4) 
and (3-14), 


/ ( 2 ) ( S ;r ) = (f - if / 2 J 1 * 1 * 1 / 2 ( 2!)" V ' 1 


(2|) s r(v + s + |i + i) 


i r a^ r (y) 


s=~ a> 


r=~co 


|)r(i * 6 -z.e) 


( 5 . 22 ) 


ro 
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/' 2 ><?; r ) = (? - e 1 '” 1 n 1/2 (2l) V 

00 _ 0° 

(21) 3 r(-v + s + n) 


^y,r^) 


=-co r (-v + 


s - r 


< (1) (E;7) - 5® 1/2 (e 2 - if /2 r‘ 1 (|) 2v (2E) v 

y«)T ‘fr 


,r, u / x 

1 V ir) 


— r f v + Lii + iVfi - E-±-i\ 

=-oo \ v 2 2/ V 2 ) 


S -v-- ! Kir) 




(2«) S (|)' 2S i S £ 


/ P (- V ~ 

r=-co 


1 \, r (7) 

f* + tX 


1 + Z-ZS. 

K 2 > 


the summation index s assumes only even values. 
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5 «7 Connection between the X- and -Z -Functions 

Due to the equality of the characteristic exponents V 
and -v - 1, respectively, in equations ( 5 * 23 ) and (5*24) or 

equations (5*22) and ( 5 * 25 ), respectively, X*^ 2 ^ and -2^"^ on 

the one hand and X^ 2 ^ and -JS^^ on the other only differ every 

-V-l V 

time hy a factor independent of | . Thus one may equ.ate 

-£ (l) (t}7) = i sin (v - n)rte" ( * i+V+l)jt ^ ^ (l) (y)^ (Sjy) (5-26) 
v it v -v-l 


^ (1) (C;r) = e (v+1)lti / (2) (7)^ (2) Ui7) 


(5-27) 


The various factors, as sin (v - and so forth, were 

introduced for convenience. "Between k ^ '(y) and k ^ (y) there 

V V 

exists the connection 


K ^ (2) (r) = - sin (v - i-Ojre”^ 1 ^ 1 (y) 

-v-1 * V 


(5-28) 


One further obtains with the aid of equations (4.l8) and (5*17) 


„h(l)/«. \ D(l)/ , -V«i 

-Z y (?;r) = (7) e COS Vrt XT 




(p+v)jti B j[ n (y + 4)^ h(2), . 

« % U,7; 


(5.29) 


COB Vlt 


# (2) (|;r) + sin w-# (l) (£;r) = e V5ti / (2) (7)x^ {2) (ur) (5-30) 


If V, u are integers, these relations are essentially 
simplified; then the X-ftmctions and the -Z-f unctions of the first 
and the second hind, respectively, are actually proportional. 
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Now the calculation of ^^(7) is left to he performed. To 

that end one may select the coefficient of any power of £ in 
equations (6.22) and (5.25) and carry out the comparison. One 
obtains 


«1 (1) (y) 


1 Vrti, v2V-2s s . . 

2® (2) ir(l + v- |i-s) 


r=-oo 


iV. (7) 

v,r 


[r( v + 3 L- + |)r(i.^)1 


r (i 

> i a (7) 
v,r 

r=-M ’ 


[ r (. v + + l)r(x . 


(5.31) 


V iV (r) 

Z. 


-2V“2s-2 

\2J FTT+T + s + 1 ) 


r=-m 


[ r (-v -4 -s + |X( 1 + 3 4- £ )J 


r=-<» 


i r a^ (7) 

v ,r 


r ( v + *** + §)4 ♦ e -t*i 


(5.32) 


Any even number is to be substituted for s in equations (5.31) 
and. (5.32). The value of ^^(7) is independent of the selected 
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special value of s. If one replaces In equation (5*32) [i by -ji 
and s by -s and then multiplies by equation (5*31), 


' k» (1) M r' 1 


(5-33) 


is originated. 


5-8 Wronski's Determinant 

Wronski's determinant of the -Z-functions of the first and 
second kind are defined by 


w„ s ^ (l> 

■Z v 


(£;y) A-z 4 ^(£;r) - -a' 

d| v V 


-J»< 2 ) 


tt;v) 4< (1) (C;r) 

d5 v 


(5-34) 


From the basic equation (2 .4g) there follows in the known 
way that Wronski's determinant of any two of its solutions is 
proportional to (| 2 - i) 1, The factor of proportionality is 
determined by substituting their asymptotic series for the 
-Z- functions of the first and second kind; it is sufficient to 
limit oneself to the first term (5-19). There results 


W 7 = - 

£J <V 


r .2 . 


Z. 


1 2 


iV 


r=-oo 


v,r 


(?) 


(5-35) 


Wronski's determinant of the X- functions of the first and 
second kind 


W, 


x v <1)(l;r) A x v t2)<s »> • ^ < 2 ) ®;r) frX^dir) (5.36) 


d| v ' v ~ ' d | 

results from Wg by using equation ( 5 . 30). 
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(5-37) 


(5.38) 


Simplifications result for the important special case p ® 0. 
First, one agrees upon omitting the index p when it has the value 

zero. Notr there is valid h (7) = a (7) and further, according 

v,r v,r 

to equation (4.1), because of P. n (l) = 1, 

00 

X (l) (l;7) = V i r a ( 7 ) (5*39) 

v v,r 


Thus, one can also write for Wronski's two determinants 


x(, 1) d;7)x( / 1) (l;7) 

% = 2 55 

I - 1 


(5-40) 


5-9 Other Series Developments of the Solutions 
of the Basic Equation 

Niven (reference 1) investigated series developments of the 
following form (the functions represented by them are called 
V- and W-f unctions) : 

> 7 ) = (w 2 + 7 2 ) 1//2 w 1 i r / r (7)t v+r (v) (5*41) 

r=- ca 3 
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w;r) » 

The relation of the variable w to I and t, is: 
w 2 -l ; 2 - y 2 . = r 2 U 2 - i) 


( 



(5-42) 


(5-43) 


For the coefficients c r and there results again a 

three-term recursion system which can he transformed into 
equation (4.4)-. If the coefficient of the principal term is set' 
equal to a Q , the equations 


l-i 

c V,r 


(r) 


r (v- + ■ y r . y , t + i) 


V,r 


(y) 



(5-44) 


(7) 

V,r 


, (v + r + ja + 1 
\ 2 




(7) 


r (y-* § - I* + i)r(i±|il) V,r 


(o) 

= 1* a n (r) 

p“(o) v ’ r 

“V ' 


(5.45) 


are valid 
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The series (5 .41) and (5*42) converge uniformly In each closed 
domain given by ( 7 ) < |w| < k>; or, expressed in the | -plane : 

1 < I'Ji 2 - l| < The bounding curve \/| 2 - l| = 1 is a 
lemniscate . Equation (5*42) is, as will be shown later, a limiting 
case of a general development, which still contains an arbitrary 
parameter and which yields as a further limiting case the 
series ( 5 . 1 ) and ( 5 * 2 ) of the -Z-functions . 

One can immediately give further series developments of 
V- and W-functions j to this end one has to replace the 
functions t v+r in equations (5*41) and (5*42) by ny +r , or the 

indices V + r by -v - r - 1, or p by -p, or one has to make 
two or three of these substitutions simultaneously. . One thus 
obtains a total of eight Y-f unctions and eight W-functions . Their 
properties will not be investigated here more closely; it should 
only be mentioned that all of them also can be expressed linearily 
by the functions of the first and second kind which is done in 
the simplest way with the aid of the asymptotic series. 

Whereas the asymptotic series of the ^-functions progress with 

powers of the asymptotic series of the V- and W-functions 

one obtains from equation (5-41), and so forth, by substitution of 
the asymptotic series of the cylindrical functions, contain powers 

of (j ; 2 - 7 ^) ^ ' , that is, (i 2 ~ l) • According to a 
suggestion by Wilson (reference 7) one can now also set up asymptotic 

series which progress with powers of (It 1) " . They have compared 
with the series ( 5 - 12 ) a slight advantage insofar as a three-term 
recursion system results for their coefficients. Correspondingly, 
for the solutions of the basic equation also developments in terms 
of cylindrical functions with the argument ? t 7 = 7(5 t 1 ) may 
be given, of the form 



(| - l) li / 2 (| + 1 ) "(V 2 t (£ t 7 ) 

t— v t v+t 

00 3 


(5-46) 


where t runs through all integers, the odd as well as the even 
ones. Y v+ ^. can again be replaced by n v+J .> and so forth. 

These developments will, however, not be followed up here. 
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6. CALCULATION OF TEE COEFFICIENTS OF THE SERIES DEVELOPMENTS 
.... IN .TERMS OF SPHERICAL AND CYLINDRICAL FUNCTIONS 
.6.1 Continued Fraction. Developments 


The solution, of the recursion system (4.4) which for r — > » 
has the behavior at infinity (4.7), can he represented by the 
convergent (reference 13) continued fraction 


a r _ 7 *3y/0 r y Pr^r+ g/^r^r+2 7 l D r +2^is-4/^r+2^r+4 

vs' 1 1 "I i ' 


(6-1) 


The solution which has the behavior (4.7) for r — > - <» can be 
represented by the convergent continued fraction 



. . . ( 6 . 2 ) 



The siibnumerators of both continued fractions are in each finite 
closed domain of y- and A, -values for sufficiently large values 
of r in the case (6.1), of -r in the case (6.2) smaller than 
one-fourth^ thus, according to a theorem on uniform convergence 
of continued fractions, the continued fractions (6.1) and (6.2), 
respectively, are in each domain of this kind for sufficiently 
large r's and -r's, respectively, uniformly convergent and 
are therewith regular analytic functions in y , and X, since the 
individual approximation fractions are functions of this kind. For 
not sufficiently large r's and -r's, respectively, then follows, 
that these continued fractions are also analytic functions which, 
however, need not in every case b.e regular « 

A solution of the recursion system (4.4) has now to be found 
which shows the behavior at infinity (4.7) for r— > oo as well as 
for r— > Then' the value of a r /a r _2 calculated from 

equation (6.1) must equal the value of this expression calculated 
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from equation (6.2). An equation results which for given V and p 
allows calculation of the separation' parameter A as a function 
of 7 . If 7 and A - v (v + 1) both are sufficiently small, the 

solution A v ( 7 ) is a regular analytic function of 7 which 
assumes for 7=0 the value v(v + 1). Thus the A^ ( 7 ) as well 
as the a r / a 0 can developed in power series in 7 with non- 

vanishing radius of convergence the magnitude of which will not he 
investigated here more closely. 


6.2 Method for Numerical Calculation of the Separation ' 
Parameter and the Development Coefficients 

The representation of the coefficient a T by continued 

fractions is also for larger values of 7 still particularly 
suitable for the numerical calculation of the separation 
parameter A and the a r 's . Mostly v, u, and 7 are given. 

Then the values p r , q r , and + A can be calculated 

numerically from equation (4.5). One starts from a value for A 
which is assumed as close as possible to the actual value and 
calculates for a selected fixed r the expression ^+ 2 /^ 

from equation (6.1) as well as from equation (6.2) . Then one 
repeats this calculation with a slightly altered value of A and 
examines whether thereby the agreement of the two values ^+2 / a r 

is improved. By further variation of A one can finally obtain 
an agreement of arbitrary accuracy. Therewith one can find the 

u, . 

value Ay( 7 ) with any desired accuracy. 

One more investigation has to be made: whether the solution 

thus found for y - 0 goes over continuously into v(v + 1 ), 

that is, into A^(0) and not perhaps into A^ (0) : for A^ (0) 

also is a solution of the present problem as can be recognized 
from the fact that equations ( 6 . 1 ) and ( 6 . 2 ) contain the values v 
and r only in the combination V + r. This question cannot be 
decided unless one has already a general picture of the 

functions A^(y) as it is given in figure 1 for n = 0 , 

v 2 1 
are integers, and real 7 s. 


v 's that 
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The number of terms of the continued fractions (6.1) and (6.2) 
to be included in the calculation corresponds to the desired 

accuracy . For large j r | the partial fractions 7 S? r < 3 r+ 2 / 0 r 0^2 

assume the order of magnitude y^/ (l 6 r 1+ ) ) thus the index r ' of 
the last partial fraction to be included will haye to be selected 
at any rate larger than j y j / 2 . 

The calculation of the a (?) 's is made by taking the 

v ,r 

value found for \^(y) as a base, and calculating a r+2 / a r from 
equation ( 6 . 1 ) and therefrom a^g/a^ a+^/a 0 , . . . 


6.3 Power Series for Separation Parameter 
and Development Coefficients 

For the numerical calculation of the separation parameter and 
the development coefficients one can for small values of | y\ make 
good use of the power series developments in terms of y . If one 
limits oneself in 'these to the first terms up to the fifth power 
of y 2 , inclusive^ one obtains, in general, still quite useful approxi- 
mations up to about j y 2 j = 5* Therefore, following the power 

series for the l^(y) shall be calculated explicitly to 7 , 

V g 

inclusive, for the to 7 , inclusive. Therewith one more 

series term is obtained than by Niven (reference l) j compared with 
Niven's cumbersome treatment, the calculation is essentially 
simplified. 

For the limiting case 7=0 there follows from the recursion 
system (4.4) 

Jx^(y) + (v + r + 1) (v + r )j a r “ 0 (6*3) 

The case where all a r disappear is not of interest since it leads 
only to identically disappearing solutions of the basic equation. 
Thus there becomes 4 ^( 0 ) = v(v + l), = 0 for r jl 0 . 
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For 7=0 all 0 r with the exception of 0 Q have nondisappearing 

limiting values. From the continued fraction (6.1) one can draw 
the conclusion • 


r - ^[l + 0(r 4 )] (r - 2, 4, 6, . 


r -2 


and from that further 



(r = 2, 4, 6, 


(6.4) 


(6.5) 


If one takes the next partial fraction into consideration as 
well, there results as the next approximation 


i r q 2 q 4 ‘ * q ;.' 

a 7 » *• * 


p,q s 


+ 7~l j~;r x " + • • • 
v..' J 2/ 1, 4 V 1 -;/ o 


p^H-r-i-S 


0 ( 78 ) 


(r = 2, 4, 6, . . .) 


( 6 . 6 ) 


Accordingly, one obtains 


a v 


-r p -2P-4 • • * 
'^2^-k * * * 

SSj * «’*> 


1 * y + \w 


9-2P-4 Q-4P-6 


P-20-4 0-4 0-6 


(r = -2,-4, -6, . . .) 


(6.7) 


One now substitutes ag/a,., from equation (6.6) and &-2/ a o 

from equation ( 6 . 7 ) into the equation r = 0 of the recursion 
system (4.4) and obtains 
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4l 


0 


P o q 2 


. , 4 p 2 q 4 

1 + 7 m 


V-2 




. * 4 q -2 p -4~ 

1 + 7 escj 


o(r 8 ) 


( 6 . 8 ) 


n . 

This equation permits the calculation of X (y) as power 

series in terms of y up to the power y^°, inclusive. At first 
one can see, "by having y approach 0, that $ 0 = O(y^) . Therewith, 
however, 0 r also is known for any r with the exception of terms 
of the second and of higher powers in y2. If one now inserts 3> r 
in this approximation on the right side of equation (6*8), 

becomes already correct up to the third power in y 2 , inclusive. 

If one repeats this procedure with the new values of the / r which 

are correct up to y^, inclusive, there results finally <f 0 and 

therewith (y) exactly up to ylO, inclusive. The performance 
v 

of this calculation as well as the calculation of the a r 's is 

not particularly difficult, therefore the results are given 
immediately . 

In order to make the representation clearer, the following 
abbreviations are introduced: 


«’r = ^rC 1 + ^ 8 r) + 

where 


( 6 . 9 ) 


D r = r (2 v + r + 1) 

5 2(4u 2 - 1) 

r (2v + 3) (2v - l)(2.v + 2r + 3) (2v + 2r - 1) 


( 6 . 10 ) 


p _ P 0 q 2 
p 2 " — > 

d 2 


Pi - V J1 4 


P-2 


. P-2 q 0 

P-2 ' 


P-4 


P-4 q -2 

P-4 


( 6 . 11 ) 
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A i “ * 2 ^ 2 * + :P -2 5 -2 i 


P 2 4 P -2 , 

B i ’ 5 2 + 5^ S -2 


0 . !£i ♦ L-£:4 

d 2 d_ 2 

1 s= 0 9 1, 2, • • • 


>- 


Then there becomes 


^(r) - v(v - i) - r 2 ^ gv 4 g 3) ^ ^ ^ - ^ ir " V 4 * V* 


(j>oB 0 - '"2 - o)r 8 + 

P -? P -h/ \ 

-( 2B -2 +8 -0 


D 


-2 


a 3 ' 2A o B l " B o A l T 3 


■ylO + o(yl2} 


P 2 P 4 


a r 


2 2 q 2 


r 


a o = 7 Dc 


4/1 e x !,4. ~_t°) 


2 J, p„ . 4f» 2 . 4 ° 

-y - r s e - 7 Vf2 + -TJ-; 

- L -n 


2A 0 8 2 + A-j_ Pi,. 


Dr 


— 


)-B 2 3 

" 




l r + o(r 10 ) 


f4 4 2s q 4 
a o“ r E 2 D 1# 


1 - 7 2 (b 2 + 6,,) 


'6 


. 4/ r 2 . - » . - 2 , x 4 . ^6 " "o’ 

+ 7 (Og + + S 4 * 


p 4 - A o + p 6 - A ~ 

Dg Dlj. 


( 6 . 12 ) 


( 25 2 + 54 ) 

(6.13) 

(6.14) 


+ 0(7 10 ) (6.15) 
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% m ^2.%% e 

a o d 2 d 4 d 6 7 


1 - 7 2 (§2 + .% + agj] + 0( 7 10) (6.16) 


a , !gw,8 a + 0(yl0) 

a o %¥s 7 } 


(6.17) 


!i sy 2-!; 


*1 _ y 2g 

JL 

B 2 
v -2 

? -4 " • 

a o " D. 

-2 

r 7 °-2 

+ 7 \ 

P -2 

j. .,6 

2A 0 I5_ , + ^ 

>d 

t 

-S' 

lOo i 


T T 


p -2 

‘ P ~2 

(2o_2 •+ 

S -4j - 




5-; 


J J 


f + o( r 1D ) (6.18) 


H„ r 4 P-2P-4 

a o D-2 d -4 


1 - 7*” (& -2 + & -4) 


+ 7 4 (s_ 2 £ + 5. 2 &.4 + 5_4 2 + A ° + -“- 6 - ^ 


D 


-2 


D 


+ o( 7 i°) (6.19) 


z£ = 6 p -2 p -4 g - 6 

a o D. 2 D_4D. 6 


1 - r 2 (s. 2 + s. 4 + a.g) 


+ 0 ( 7 10 ) ( 6 . 20 ) 


1* . 7 s !rgl.€-6 p -8 + ( 10, 

<o B. 2 D. 4 E.gD. 8 ) 


( 6 . 21 ) 


For the case excluded, above where v has fractional values of 
one-half, the convergence radii of these series equal zero. It seems 
therefore probable that the convergence radii are functions of v 
which can be infinitely large for special cases, but not in general. 
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6.4 Power Series Developments 

Since occasionally power series developments of the solutions 
of the basic equation (2.4g) also can be useful, they will be 
briefly discussed below. 

One can of course obtain them at once by substituting in the 
series (4.1) and (4.2) for the X-functions of the first and second 
kind, the known power series developments of the spherical functions 
in terms of powers of I; one thus obtains power series for the 
solutions of the basic equation which converge in the circle }|( <1. 
The problem of the Laurent- series for 1 < III < “ need not be 
discussed further since they are already calculated in. equations (5*22) 
to (5.25). However, one can obtain these developments directly. 
Therewith a new method for the calculation of these functions and 

particularly of (7) is found. 

One starts from the differential equation (2.10) which is 
written in terms of | rather than of £ . 


<» e - a ♦ 0 S 


+ r 2 S e - 

I 


v = 0 (6.22) 


For the integration one tries the statement 


, % .a fc v+s 
8 (r)i I 
v,s 

s=-°° 


(6.23) 


Then there results for the g 8 (the indices p and v as well as 

the argument y in general are again omitted) the three -term 
recursion system 
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(v + s - fi + 2) (v + a - |i ■ 

+ [(v +3+1) (v *>• s) “ 

(s a 0, ± 2 , t4. 

There exists a solution with the 

for s— >°° and a solution with the behavior at infinity — — — > 1 

e s-2 

for s — » - oo . The quotient of the two solutions is independent of s 
only then when X assumes certain distinct values . As one can see 

by comparing with equation (5.23), these are just the values X^(y) . 

From the behavior at infinity of the coefficients g one can 

s 

conclude at once that the series (6.23) converges in the 
domain 1 < 1 1 1 < 

If one substitutes in equation (6.24) for the coefficients g B 

the coefficients calculated already in equation (5.23)/ there results 
after elementary transformations 


ijgg+r, 

m0 

. . .) (6.24) 

behavior at infinity > ^ 

So a 


* r <,r (y) 


(2y + r 


+ s 


r=-» r (' v - t i^ + i ) r ( 2 + E i^) 

1) (2 + r - s) Jjv + r + 1) (v + r) - xj 

*1 

- 7 2 (v + 3 “ |i - 1) (v + s - |i) 

(s = 0, +2,' +4, . . .) 


■}- 


(6.25) 


These relations can be used, like equation (5-15 ) , for the control of 

u , 

numerically calculated values of the a (7) . 

v,r 
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The recursion system, (6. 24). is, except for the case of v,u 
being integers with v ^ | pj s? 0, probably less suitable for the 
numerical calculation of the 1^(7) than the. continued fractions (6.1) 

and (6.2) . 

Ordinary power series with increasing powers of I result 
for the solutions of the basic equation if one sets egual 
g-2 = s -4 “ * * • = 0 an<a - feq'uires = 0. 

Then there results for v the determining equation 

(v - |i)(v - |i ” 1) = 0 (6.26) 


Therefore ■ v has here a meaning different from the one it had so 
far. The behavior at Infinity of the g s for s — - » is simple: 

all of them disappear. The behavior at infinity for s — > <» is 

^ . . 

given by — — >1 or — $> The first case is the standard 

6 s-?_ s 

case: the power series converges for f|| < 1. The second case is, 

for v and p being integers with' v = fpf = 0, realized for a 
solution of the basic equation, the X-function of the first kind; 
the power series then converges for all finite I . 

It will be best to make the numerical calculation of the 
coefficients of these power series which are convergent in the 

JJL 

unit circle so that first \ (7) will be determined according 

v 

to the method given in section 6.1, or, for smaller values of 7, 
from the series (6.13); the coefficients g 8 can then be calculated 

from equation (6.24) for- each of the two v -values given by 
equation (6.26). A special but simple problem will then be left: 
how the two calculated power series are connected with the X-functions 
of the first and second kind. 
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7. EIGENFUNCTIONS OE THE BASIC EQUATION 
7.1 Limitation to v,p Being Integers ; v = |p| = 0 

The determining factors for the eigenvalues of the separation 
parameters \ and p and, if occasion arises, of the wave coef- 
ficient lc, are the domain of space which was taken as a "basis 
and the boundary conditions on its "boundary. This treatise is 
limited to the most important type of eigenvalue problems of thiB 
kind; for them the domain of space lies either within an ellipsoid 
of revolution, or between two conf ocal ellipsoids of revolution or 
outside of an ellipsoid of revolution. The first two cases will be 
called problems of inside space, the last case problem of outside 

space. The entire domain -1 = tj =5 1, 0 I 2 cp Is 2jt becomes then 

effective for the two coordinates t) and cp. Boundary conditions 
in t| and cp do not appear then; they are replaced by the require- 
ment that the wave function for T) = tl remains finite and that it 
is single valued, that is, that it has the same value for cp + 2jt 
that it has for cp. The latter requirement leads to p's that are 

integers, the first one to v’s that are integers v = |p| - 0. 

That the X-functions of the first kind remain finite at the points 
q = tl follows directly from the series (4.1) by taking the 

estimate i (t ) = * | £ + \fT^ - l| , which is valid for 

I V v , 1 

this case, as a basis. 

Following, n will always be written for v and m for p 
where v and p are real integers; for the present, n ^ m 0 is 
assumed. The case of negative m's, the absolute amount of which 
is = n, iB then obtained at once from equations (4.15), (4.17), 
and (5*20) . 

The calculation of these special functions was practically 
settled amongst other things in the last sections; even though it 
was assumed there that neither v + p nor v - p are integers, 
almost all results can nevertheless be taken over as simple limiting 
processes demonstrate. Only a few particularities result, compared 
with the general case; they will be discussed below. 
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7.2 Breaking Off of tho Series 

If q , = 0 for a positive r' or P r i = 0 for a negative r', 
the a^'s break off to the right or to the left, that is, 


for q , = 0: r'>0: a , = a , „ = a 


V 


for p , 
x r 


= 0; r' < 0: 


r 


a t 
r 


r'+2 r'+4 


— a 1 ~ — a 


V-2 


'-4 




= 0 


= 0 


(7-1) 


is valid which follows in the simplest way from the continued 
fraction developments (6.1) and (6.2). These .cases occur when p - v 
is a positive integer or when p + v is not a negative (Sic!) integer. 
Since it was presumed 0 = m = n the first possibility does not occur, 
but the second one does always occur, that is, for all admissible m,n. 
Here again two cases must be distinguished which are both originated 
from p r i = 0: 


n + m+ 2=5-r’>0 or n + m+l = -r'>0 (7 *2) 


In the first case m t n is an even number, in the second, an odd 
n umb er; a r >+2 ~ s first nonvanishing a^. for the b^ ' s there 

follows from equation (4.12) that, they disappear for all r = r' + 2m. 

Further, all P m (| ) disappear for n + r = -m, -m + 1, . . . m - 1. 

n+r 

The developments of the X- and %-f unctions of the first kind begin, 
therefore, for n - m = even with 


-^ a . n+2 

< ( 1 ) (Ur) - (f - yf'\- 


.m' 


'"^(5) + • • • 




r ( 7.3) 
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and for n - m = odd with 


m(l) . m m-n+1 m . . 

X n <«»>..- a n,m-n+l 1 WO 


+ a S,m-n + 3 im ‘ n+3 ^ (l)+ • • • 


■mi 3 


- (f ■ fY'* i WU 


> ( 7 - 4 ) 


+ b m f (t ) + 

n,m-n+3 y Btf3' b 


The series for converges for all finite The 

corresponding formulas for the ^-functions of the second kind result 
if the functions t v+r (£) are replaced "by n v+r (£)* The developments 

for the X-f unctions of the second kind show a special "behavior . The 
spherical functions of the second kind "belonging to the vanishing 
coefficients a^i, a r i_ 2, • • . "become infinitely large in such a 

manner that their products have finite limiting values . The 

coefficients a m (7) are defined "by 
n ,r* ' J 


lira 



m _m 
a n,r ~n+r 


(I) 


for m = 0, 1, 2 } . . . and v + r + m — > - 1 , -2, . . . 


( 7 . 5 ) 


Then there becomes 


m . . m+n+1 , m 

a n,r'" ~ lim T (1 + v + m)r(- v-m)a^ ^(7) 

for m = 0, 1, 2, . * . and v + r + m — > -1, -2, 


and v + r + m 


( 7 . 6 ) 
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and the series (4.2) reads 


£ (a) (Ur> -J_ 


n 


.r m 
i a 


«*£,»>♦ 2 _ 

r=r '+2 


n,r 


r=-os 


( 7 - 7 ) 


The series (5.41) for (w ;y) breaks off only when n - m 

is an even number; for odd values of n - m the coefficients c^ Ay) 

have the indefinite value oo.O if r = r 1 with finite limiting 

value. The series (5.42) for (w;y) , on the other hand, breaks 

off only when n - m is an odd number; for even values of n - m 

the coefficients d m (y) have the indefinite value » .0 with 

n,r 

finite limiting value if r = r . Similar conditions exist for 
the other V- and W-f unctions . 


7.3 A Few Special Function Values 

From the series (5>l) one obtains when arg (£ 2 - 1) = n 
for £ = 0 ■ ' 


m(l). 

* (0 >7) = 

n 


lift i m (r/ 2 ) m b^ m rn (y)j r 0 + m) for 


r 


agg ( 1 ) (o;y) 

d? 


0 


for 

n - m 

even 






for 

n - m 

odd 







for n 

- m even] 


i m (y/ 2 ) m ^n 3 m-n+l^j V (2 + m ) for n - m odd 


T ( 7 - 9 ) 
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The X- and -2-functions of the first kind are for the 
index values n and m considered here either even or 
odd functions of I or £, respectively, according to 
■whether n - m is an even or odd value. Furthermore, because 

of j^Kcos 0 + i.O) * (-iJ^Q^cos 0) - ~-P m (cos 0 + i.O) there is 
n n & * n 


y% ( 2 ) (0}r) =-|«i 2® (l) (o;r) 

<< ( 2 ) (o,- 7 ) ! ax^ (l) (o; r ) 

1 rti 


for n - m even 


for n - m odd 


(7.10) 


(7.11) 


From Wronski 1 s determinant (5*35) follows 


m(\) 3£n (2)(0;r) 

ra m(l) ( o. r) n . 


I 

rs- oo 




for n - m even (7*12) 


(o\ d£“ (l) (0; 7 ) 

yg m ( 2 )(0;r) — B - 


% 


i r b m (r) 


n,r 


r«~oo 


for n - m odd (7*13) 


Therefrom the -2-function of the second kind and its derivative 
with reBpect to £ for £ » o can be calculated at once . 


7 .4 Connection between the X- and -Z -Functions 

If v,u are integers, considerable simplifications occur in 
the relations (4.15) to (4 .2o), (4.24), (4.25), (5*6) to (5*9), 
(5.17), and (5.18). They are so obvious that they need not be 
discussed further. Equations (5.29) and (5*30) now assume the 
simple form 
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■e” (1) (5;r) - *“ (i, W^ (i) ({5r) <* - 2 ) (T-W 

For the n m ^(y) simpler expressions can "be obtained if s 

XX 

in equations (5. 31) and (5*32) is selected in a suitable manner. 

The same expressions, however, result in an even simpler way if 
one substitutes in equation (7 «l4) and in the derivative of this 
equation with respect to I, respectively, the special value | = 0. 

If one expresses ^(Ojy) and dE^ (0 }7) /&&> respectively, 

according to equations (T -8) , (7*9), (7*12), (7*13), using 
equation (5*13), there originates for n - m = even . 


m(l) . 1 l/2 ^n,m-n^ y ) 

k (7) “ 5« 1 V?J /■? \ 

n ' 1 (o;7)r(f + ^ 


y (7.15V 


m(2). . 

"n (r) ' 




.-1 C (1> ( °;r)r(I - .) 

(7) 


n,-n-m 


and for n - m « odd 


«™ (1> (r) = 1 m (|f 1 


m , . 

1 . -1 (7) 

n .m-n+l 


dx“(l)(o;7)/d| r^| + m) 


> (7.16) 


K m(2) (7) = 

n v/ ' 


-it 


•l/2 -m/^\m-2 


ax" m( l) (0;y)j dg r(|- - m) 


w 


m / % 

a n,-n-m!l ^ 
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i 


By 2^(1) ( 0;r ) 

and ax“ (1) (Ojr)/ai the values of these 

functions are understood which result when § goes towards zero 
from the positive imaginary half plane. 


The distinction between even end odd n - m can be' avoided 

/j \ 

if one sets, for instance, s equal zero in calculating the » v 1 (y) 

from equations ( 5 * 31 ) and ( 5 * 32 ) j the formulas' ( 7 * 15 ) and ( 7 * 16 ), 
on the other hand, have the advantage of greater simplicity. 


7 *5 Normalization and Properties of Orthogonality of the 
X-Functions of the First Kind 

The eigenvalues of the basic equation 1 ^( 7 ) are always real. 
Proof of it is given in the known manner. Equally simply it can be 
shown that the functions are orthogonal to each other, 
that is, 





0 


(7.17) 


is valid for n ^ n * . 

By inserting the series (4.1) into (7*17) one can also express 
this property of orthogonality for even differences n - n' thus: 



1 n,r (, ' )1 ’n' < r+n-n' (r) ST"Ir + 1 * 0 for n * 


rc- co 


For the normalization integral one obtains 


s S <1) 6 ”K (1) ® ») « - (7.19) 

"1 
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7 .6 Generalization of 3?. E. Neumann's Integral Relation 

In the case m = 0 one obtains a second solution of the basic 
equation ( 2 . 4 g) which i-s independent of Xp) (|;?) in the form of 

the integral 

F(?) f xC 1 ) (t;r) dt ( 7 * 2 o) 

2 J-i I-* n 


The fact that this integral actually represents a solution 
of equation (2 .4g) is confirmed by substitution. The calculation 
is reproduced in detail in Bouwkamp (reference 10) . For large £ , 
t and I in the denominator of the integrand cancel in first 
approximation and one can see then at once that <T(i) is proportional 
to the -E-function of the fourth kind. The integral over t can then 
be evaluated according to equations (8.26) and there originates, 
because of equation ( 5 * 12 ), 


J(!) 


-ir 


P^gJS) - i2 ( p(n;s) 


z 




(1) 


n 



00 



(7-21) 


According to equations (7*8) and (5*33) the -Z-f unctions are now 
converted to X-functions. Because of equation (5*39) there results 
finally 


x ( p(i;r) - i7 


(r) 


x^flir) = | ; 


n 


r * 1 iy(t-i) ,-v 

■ X P (t;7) dt (7 ,22) 

o-i 


Therefrom results for 7=0 F. E. Nexmiann's integral relation 
between spherical functions of the first and second kind. 
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7 .7 Zeros of the Eigenfunctions 

For m > 0 the zeros of the 'basic equation are situated 
at £ = tl, respectively, since they there have the 

behavior (l 2 - x) m / 2 . if one divides the eigenfunctions by this 
expression, the quotient does not have zeros at I p tl. In order 
to understand this, one need only enter the basic equation (2.4g) 

•with the expression (I t l) m / 2 multplied by a power series 

in (| t 1). The zeros of - l) jy) are all simple; 

if they were not simple, all higher derivatives of the eigenfunctions 
would have to disappear there also. Since it is, however, a non- 
identically vanishing analytic function, this case can never occur. 
Further properties of the zeros of the eigenfunctions follow from a 
simple consideration of continuity: namely, that in a nonsingular 

point of the basic equation a zero cannot be newly originated for a 
change of . . . * and an already existing one cannot vanish. There- 
with the problem of the number of the zeros is essentially reduced 
to the problem of the number of zeros of Legendre's and their 
associated polynomials and of Bessel's functions with an index of a 
fractional value of one - half . 


One deals first with the X^^(t;y) with real £ and y, 

that is, with the eigenfunctions of the prolate ellipsoid of 
revolution. All zeros are real; for this is valid for y = 0. 

If, namely, for a change of y a complex zero would originate, the 
conjugate -complex would originate along with it; but it contradicts 
the simplicity of the zeros, that a real zero splits into two 
complex zeros . The number of the zeros in the interval -1 < I < + 1 


equals n - m, that is, the number of zeros of P m (£) in this 
interval. The zeros outside of this interval go over Into the 
zeros of / n (y£) for y — > 0; the asymptotic distribution of 


J n+l/2' 

the zeros for large £ is the same as the distribution of J 
for arbitrary y . 


n+l/2 


(76) 


For the eigenfunctions of the oblate ellipsoid of revo- 
lution X^(!;iy) with real I and y also n - m zeros are 
situated in the interval -1 < £ < 1; but now the remaining zeros 


Translator's note: 


. . . missing in the original. 




\J1 

CT\ 


lie on the positive and negative imaginary semiaxis of the complex |-plane; they agree 
asymptotically as veil as for y = 0 with those of J n+ ^/. 2 (i 7 ^) . 


7*8 .Integral Equations for the Eigenfunctions 

The most important integral equation for the X-functions of the first kind is 
derived from the known relation 


n+m 

Qi 


(n + m) 1 
(n - m) ! 


-m t \ 

2 'M 2 ) 



u 


-i 




(7*23) 


If one writes y£ for z, replaces n by n + r, multiplies by I r a^ ^(y) 
forms the sum over all even r, there results 


and 




-m 

7 


_m(l) 


(ri;r) 


Pi 

i 


irln 

e 




( 7 * 24 ) 


Equation (7 .24) Is a homogeneous linear integral equation of the Eredholm type 
for the X- and -Z-f unctions of the first kind, respectively, with the symmetric 

kernel (l - r) 2 )(l - | 2 ) . If one replaces in accordance with equation (7 .14) the 

•Z- function of the first kind by the X-function of the first kind, one can see that 

the connecting coefficient ^ ( 7 ) in this integral equation plays essentially 

the role of the eigenvalue parameter. 
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The integral equation (7*24) can he generalized. By selecting 
another path of integration one can, for instance/ also express the 
Z- function of the second kind (that is, also the X-function of the 
second kind) "by an integral over the X-function of the first kind; 
furthermore, equation (7*24) can he generalized to the case of 
arbitrary v,u,y. However, the respective results shall not be 
discussed here. Kotani (reference 8) indicated a general principle 
for obtaining more general integral equations for the X- and 
Z-functions, respectively. 

Integral representations for the X- and -B-f unctions have not 
become known so far. It seems that the integral equations of the 
type (7 *24) or of another kind also can be substituted for them and 
replace them; thus equation (7*24), for instance, represents a 
very useful starting point for the investigation of the I -asymptotics 
of the X-f unctions of the first kind. The integral equation (7*24) 
can perhaps also be applied when the values of the X-function of the 
first kind are known only in the interval -1 < 5 < 1 and are to be 
calculated for arbitrary real and complex 5 . (Compare the discussion 
on the y -asymptotics of the eigenfunctions in the following section.) 
As Mbglich (reference 4) has shown, the integral equation (7 *24) can 
also be used for obtaining developments of the X-functions of the 
first kind in terms of powers of y . 


8. ASYMPTOTICS OF THE EIGENVALUES AND EIGEMFUNCTIONS 
8.1 Asymptotic Behavior of the Eigenvalues 
and Eigenfunctions for Large v 


The continued fractions (6.1) and (6.2) do not only have the 
property to yield a development of X.^(y) in terms of powers of y 

but in addition one can obtain from them a development in terms of 
powers of v"-*-. It is more favorable to set up a development in 

terms of powers of (2v + 1)"^, because then the odd powers 
of (2 v + 1) are eliminated because of equation (4.9). The cal- 
culation itself is relatively simple so that the result can be 
given immediately 
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x£(r) » v(v + i) + |r 2 + — ^ [t 1 * - i6p 2 ) 

v 2 • 8(2v v 1) 2L 


2 ' 

7+7 


4 


2(2 V + l) 4 


(4 - l 6 |i 2 ) 7 2 + Qji - + 0 (2V + l)" 6 j (8.1) 


Presumably this series is not convergent hut has asymptotic 
character. 

In order to form a Judgment on the usefulness of the series (8-1) 
for numerical purposes, one gives for several cases the numerical 

value of the remainder term denoted, by 0 [(2v + l) in comparison 
with the value of the separation parameter X itself . 


n 

2 

4 

6 

8 


11-7904 

25.2513 

47.10953 

77 -06246 

Remainder term 

-0.134 

0 .0132 

0.00095 

0 .00017 


In a similar way one obtains the following expressions for the 
development coefficients of the eigenfunctions for large values 
of v. 


^2 _ 
®o : 

e 

8(2y + 1) p 

Hi_ 

: (2v + 

I ) 1 

2 + 0 

(2v + 

x,-! 

(8.2) 

a -2 

a 0 

7 2 

8(2V + 1) 

p 

^7 

2(2.v + 

0 + ° 
l ) 2 

(2v 4 

• l>-3] 

( 8 . 3 ) 


a 4 7 ,+ 

0 + 

o| 

( 2 V + 

1,-3 j 

*'1 



a 0 " 32 (2 v + 

l ) 2 











- (8.4) 


a -* 

_ . . ^ 

0 

|~( 2 v + 

i)~7 




a o ' 32 (2 V + 

l ) 2 


mJ 
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8.2 Asymptotic Behavior of the Eigenvalues for Large Beal y 

One limits oneself here to n and m -which are integers 

(n = m = 0) and to real large y which may he assumed to he positive 
■without essential restriction. Thus one obtains the asymptotics of 
the eigenvalues and eigenfunctions for the wave equation in the 
coordinates of the prolate ellipsoid of revolution. 

An approximate picture of the eigenvalues and eigenfunctions 
is obtained if one puts the basic equation into the Liouville 
standard form 



cos^ 9 



sin 2 9 



- 0 


q =F 1 U)(l - | 2 ) 1/4 


I = cos 0 


> ( 8 . 5 ) 


and interprets it as a Schrodinger wave equation in the 
interval 0 ^ 6 ^ n of the potential energy (in suitable units) 

p 2 - ~ 

7 2 cos 2 0 + - — - - p (8 .6) 

sin 2 0 k 


It has for large y at 0 = | a very narrow minimum and can there 

very well be approximated by a parable. Then, however, there results 
just Schrodinger' s wave equation of the harmonic oscillator, for 
which eigenvalues and eigenfunctions are known. 

In order to obtain also higher approximations it suggests itself 
to attempt a similar formulation as in equations (4.1) and (5.1). 

One sets 


(« 8 - O " 72 


y=- oo 


(8.7) 


L 


6o 
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equal, where the D n are Hermite 's orthogonal functions and the 

functions of the parabolic cylinder, respectively. By substitution 
of equation (8. 7) into equation ( 2 . 4 g) there results, if one utilizes 
also the recursion formulas and the differential equation of Hermite 's 
orthogonal functions (see Magnus and Obcrhet linger (reference ll)), . 
the five-term recursion system . . 


d > - 4[1 d 0 + 
r -4 H r-2 


- 4 \ + 87(11 + r + =0 


hra. c 


- 2 (N + r) 2 - 2 (N + r) - 3 $ r + 4 u(N + r + 2) (N + r + l)d 
+ (N + r + 4 ) (K + r + 3 ) (N + r * 2) (n + r + l)d r+]+ = 0 


(r=0, 12, 14 , . . . ) 


( 8 . 8 ) 


The series (8.7) is probably not convergent; it rather represents 

^r 

an asymptotic development in the sense that limits — = 0 for all 
even r f 0 , or, as one concludes from that and fro3u equation (8.8) 


= 0 ( -i-r/2) 
l9 o 

Q( -l-r/2 } 
*0 

(r = 0, 2, 4 , . . . ) 


( 8 - 9 ) 


By a method of successive approximation the d r and \ can be 

represented as power series in 7 -1 . The calculation is elementary; 
thus only the result is given. It is 



MCA TM No. 1224 


61 


x“(y) » (2N + 1)7 r |(2H 2 + 2N + 3 - 4m 2 ) 

- j|^(2N + 1 ) (h 2 + N + 3 " 8m 2 ) 

+ _JL_. ^48m 2 (2N 2 + 2N + l) 


- 5(N i! ' + 2 n 3 + 8# + 7N + 3)] + 0(y"3) (8.10) 


The connection between N ; n, and m is given by counting the 
zeros. For 7 — > » the development (8.8) is reduced to the principal 
term with the N real zeros of EerrrJto’s I'T^ polynomial, whereas 
the X-function with the indices n ani m to be approximated has 
exactly n - m real zeros in the internal -1 < I < 1. Therefrom 
follows N = n - m. For negative m one inserts instead N «* n + m. 


8.3 Asymptotic Behavior of the Eigenfunctions for Large Real 7 

The asymptotic representation of the eigenfunctions results by 
calculation of the coefficients -3 v • They read, aside from the terms 

M O 

of the order 7 , 
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- ■ 357^ - 2 5H - 360 ' 


4 -2 

_ JL_ 

*0 ~' r 


*4 

*0 

^ -4 
*0 


1 + ~1_(N 2 + 27 N - 10) N(N - 1) 
327 


_ - 1 

1 4 - JL/n + 2 - 2m 2 ') 

327 

27V 2 JJ 


> ( 8 . 11 ) 


_ 1 

1 + 

- a + sJ) 

" 327 

27 V 

2 )\ 


N(N - 1) (N - 2) <N - 3) 


d 


-6 


m 


m N« 


* 0 1287 £ 


’8 


■3 0 1287' 

fl 


2 '(N - 6)1 


-8 


A 0 20487 ' 


2 ' 


Q 20487 


1 _ N) 

2 Ts - 8 ) ! 


According to the type of derivation, however, the eigenfunctions 
are approximated "by these series only in the interval -1 < g < 1 . 

In order to obtain an asymptotic series also for other | one starts 
from the integral equation (7*24) and substitutes for the X-functions 
in the integrand the series ( 8 . 8 ). Therewith the asymptotic develop- 
ment of the eigenfunctions for all g is knownj in particular, their 
behavior can be investigated where, besides 7 , I also is very 
large . Since now the eigenfunction for all g is asymptotically 
known, one obtains the solutions of the second kind by calculating the 
integral in equation ( 7 .24) with the asymptotic series of the eigen- 
function and by means of another appropriate path of integration. 
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The zeros of the eigenfunctione located in the interval -1 < I < 1 



crowd for large 7 more and more around 5 * 0; 


this, one has only to divide the zeros of Hermite ' s 
hy \f£y and therewith to convert to I -scale. 


in order to understand 
polynomial 


The domain of validity for equation (8 .10) and (8.11) extends over 
the indicated domain; thus originates, for Instance, for m = from 

equation (8 .10) the asymptotic representation of the eigenvalues of 
Mathieu's differential equation found by Ince (reference 14) . However, 
the limits for this domain of validity shall not "be submitted to closer 
investigation here . 


8.4 Asymptotic Behavior of the Eigenvalues for 
Large Purely Imaginary 7 

One limits oneself again to n and m that are integers 

(n = m = 0) and to purely imaginary 7 of large absolute value. 

This procedure yields the asymptotics of the eigenvalues and eigen- 
functions for the coordinates of the oblate ellipsoid of revolution 
(reference 10) and for the so-called inner equation for the separation 
of the wave equation of the ion of the hydrogen molecule (reference 7) • 

The method applied in equation (8.2) fails here; namely 

would become purely imaginary and the D n (Z\f 2 y ^ would, for large £, 

no longer decrease exponentially, but increase exponentially; they 
would, therefore, be no longer appropriate for the development of 
the eigenfunctions . The wave mechanical picture of the differential 
equation (8.5) shows that in the case of purely imaginary 7 two 
domains with low potential energy are present at 0=0 and 0 = 2jt, 

which are separated by a high potential peak with the maximum at 0 = 

One may, therefore, expect beforehand that the eigenvalues will 
degenerate in first approximation; their split-up is exponentially 
small in 1 7 J ; it is the larger , the higher the eigenvalue . For each 

eigenvalue there is an eigenfunction symmetric with respect to 0 = 3., 
that is, |=0 and an asymmetric eigenfunction. 

The mathematical treatment is as follows. A singularity is made 
to move to infinity. Then one obtains from equation (2.4g), aside 
from an elementary transformation, ..the differential equation of 


L 


Klcu 
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Laguerre T s orthogonal polynomials. This suggests for the solution 
of equation (2.4g) the formulation of Svartholm (reference 7) 


*1.(8) * (l 


t=- 03 


(m) 

:+t UPC 1 


-D] 


( 8 . 12 ) 


wherein i7 was set equal to p; again it does not mean an essential 
restriction if p > 0 is assumed. By substitution of equation (8.12) ■ 
into the differential equation 2 , 4g) , application of Laguerre 1 s 
differential equation, and the recursion formulas for Laguerre’ s 
polynomials (compare Magnus and Oberhettinger (reference 11)), -there 
originates in the known way for the ct a three— term recursion 
system. With the abbreviations 


X t= -p 2 + 2rp - ~(t 2 + 1 - m 2 ) +.A (8.13) 

, ' t = 2N + m’+ 1 ; (8.14) 

4A t = (t + 2t - l) 2 - m 2 ; = 2t(r - 2p + t) (8.15) 


the recursion system reads 


°t+l A t4l + ff fl A t 


(A + Pt) a t ^ 


“Nj -N +1^ -N + 2j 


(8.16) 


Therefrom follows for A the transcendent equation 


A = 


All 


A 2 


A, 


2 i 


A.i 


1 - a — i . . . + 1 - 

A + P^ |A + Pg A + P.. ^ A 4* 


(8.17) 


from which A can be obtained as series in terms of powers of p“^. 
Therefrom then results 
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X™(7) “ ‘P 2 + 2rp " ~(t 2 + 1 - m 2 ) - --(t 2 + 1 - m 2 ) 

11 ^ Op 

- ~2 jpr4 + 10r 2 + 1 - Sb^Ot 2 + 1) + ffl^l 

- + H^t 2 + 37- 2h^(23t 2 + 25) + 13^ + o(| 7 i -4 ) (8 .18) 


8 .5 Asymptotic Behavior of the Eigenfunctions 
for Large Purely Imaginary 7 
For the coefficients of the development (8.12) 

+ 1,2 ' ^ - Hr) + Hr 3 ) 

sj = - 1)2 • m2 ] [< t - 3 > 2 - - °(m " 3 ) 

if * [ (t * 1)2 ■ m2 ] C 1 - H 1 ) + °(m 3 ) 
sf ‘ 5i?[ (T • 1)2 • mS ] [ (T ■ 3)2 • 4 * “(w 3 ) 


(8.19) 


is valid. 

The significance of t and N results again from counting the 

zeros. The principal term L^ [-217(1 - ,|)J has N zeros which 

for large (7^ lie all closely to gal. The real eigenfunction 
has again N zeros in the neighborhood of g a -1. For odd n - m 
another zero at I a 0 is added. The sum total of the zeros n - m 
equals, therefore, 2N for even n - m and 2N + 1 for odd n - m: 
thus 
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t =» n + 1 = 2N .+ m + 1 for n - m = even 
t = n = 2 N+m+l for n - m = odd 



(8 . 20 ) 


is valid. 

Baber and Hasse 7 (reference 7 ) calculated the series (8.l8) -with 
the exception of the last two terms; only for the special case N = 0 
they give also the last two terms; Bouwkamp (reference 10) calculated 
the series (8.l8) with the exception of the last term for the special 
case m = 0 .- The asymptotic series (8.l8) can still he used 

for m <= t~f it then goes over, exactly like equation (8.10), into 

the asymptotic series for the eigenvalues of Mathieu's functions 
(reference Ik) . 

For large values of \y | the eigenvalues move closer and closer 
together in pairs so that the asymptotic series (8.l8) for the eigen- 
values of each pair are the same (see equations (8.20)); that is, 
the difference of the two eigenvalues has a stronger tendency to 
vanish with increasing 1 7} than any power of l/J 7 | . 

(Compare table 11 .) 

The series ( 8 . 12 ) for the eigenfunction is useless in the 
interval - 1^1 =■ 0 . There an approximation must be attempted 
starting from the point | = - 1 . Since the eigenfunctions become 
exponentially small in the neighborhood of g = 0, one can build 
up the eigenfunction in the entire interval -1 < | < 1 by combi- 
nation of the two approximations starting from -1 and 1 and one 
obtains 


x“ (1) (Si7) = Constant (l - | 2 ) m/2 I^tM 1 ' 6 ! 

t=-00 L 

1 °' ps 


(8.21) 


For even n - m the positive, for odd n - m the negative sign is 
to be selected; in the one case the eigenfunction is symmetric, in 
the other antisymmetric with respect to the point g = 0. 
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'*What was said in section 8.3 is valid for the asymptotic 
calculation of the eigenfunctions and the functions of the second 
kind for any complex | as well as for the limits of the domain 
of validity in the variables V,n,y of the asymptotic representations 

In order to show the us© of the asymptotic series for numerical 
purposes one compares for m ■** 0 a few eigenvalues with the values 
resulting from equations ( 8 . 10 ) and ( 8 . 18 ) by giving the value of the 
remainder term 0(y"3) and 0 (p*’^'), respectively. 


n 

. q 

r ..1 J.utpi . .. 
i \ 

0 

0 

2 

7 2 

10 

-25 

-100 

-100 

X n ( 7 ) 

2.305 

-16.07904 

-81.02794 

- 45 .48967 

Remainder term 

-0 .025 

-0.0l6l6 

-0.00008 

-0 .01528 


9 • eigenfunctions oe the wave equation in rotationally 

SYMMETRICAL ELLIPTIC COORDINATES 

9.1 Lame's Wave Functions of the Prolate Ellipsoid of Revolution 


By separation o.f the wave equation in the coordinates of the 
prolate ellipsoid of revolution one obtains the following solutions 
of the wave equation 


u 



( 9 . 1 ) 


A, B, C, D, E, F are arbitrary constants, v end u arbitrary 
real or complex parameters; the significance of 7 is given by 
equation (2 . 8 ) , thus 7 is real . The coordinates | and £ *= 7 g , 

respectively, n and cp are real as well. Under 3* x ^ 1 » 2 ^ (ij‘j 7 ) one 
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understands (r| + i x 0;r) . According to the kind of the 

boundary value problem presented, the arbitrariness concerning the 
constants and parameters is limited; then such solutions of the wave 
function have to be determined which remain finite for the entire 
domain of the eigenfunctions. 

Following, as before, three-dimensional domains only are dealt 
with which lie inside or outside of an ellipsoid of revolution or 
between two confocal ellipsoids of revolution. Then the domain of 
the coordinates i} and is given by -1 = 1=1, 0 = cp = 2jt. 

The requirement of single -valuedness and finiteness of the eigen- 
functions then leads to v = n, n = m, n *= |mj 0, and D = 0 . 

The eigenfunctions are written in the form 

(n = 0, 1, 2, . . . ; m - 0, tl, +2, . . ., In) (9.2) 

The domain of variables in i is denoted by ip a I Is ig 

and £p - £ - £g, respectively . For the prolate ellipsoid of 

revolution there is always 1 ts |^. For inside space 

problems ig = finite, for outside space problems infinite. For 

inside space problems boundary conditions for ip and |g are to 

be prescribed. This results in two linear homogeneous determining' 
equations for A and B ; they can be satisfied only for certain 
distinct values of 7, that is, for certain eigenfrequencies; in 
that case they fix the ratio A:B. In case ip = 1 a boundary 

condition can be prescribed only for ig > 1; the boundary condition 
for ip = 1 is then replaced by the requirement of finiteness of the 
eigenfunction at the singular point ip =1; it leads to B = 0. 

For outside space problems the boundary condition for Ig = » 
is eliminated; the functions (9*2) have for i 2 — >00 for arbitrary A 

and B an oscillating behavior. One can see that immediately from 
the asymptotic series (5 .12), The boundary condition at i = ip gives 

the ratio A:B. For ip = 1 this boundary condition in turn is 


u“(M,<Pjk) = 
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eliminated and B ‘becomes B = 0.. A condition for the frequence 
does hot exist j all wave coefficients are admissible, the spectrum 
is continuous and extends from k = 0 to k = 

9 .2 Lame ' s Wave Functions of the Oblate Ellipsoid of Revolution 

The solutions of the vave equation originating by separation of 
the vave equation in the coordinates of the oblate ellipsoid of 
revolution are obtained from equation (9*1), by replacing y there 
by ti y. Here also only the three-dimensional domains characterized 
in section 9 *1 are dealt vith and the eigenfunctions can, therefore, 
be -written in the form 


u^(i,q,cp;k) 


n 


(n « 0, 1, 2, 


Azg (l) (£;i y) + Bz“(£;i r) x“ (l) (ii;i 7 )e imcp 

. . ; m » 0, ±1, ±2, . . . ,±n) (9-3) 


The domain of variables in q and cp is the same as in the coordinates 
of the prolate ellipsoid of revolution. The domain of variables in £ 

is again denoted by £-^ = £ = £2 • ‘ For the oblate ellipsoid of 
revolution there is 0 = £ 1 . What was said in section 9*1 
for £p>l is valid for inside and outside space problems 
vith £1 > 0. However, whereas there tp “ 1 was a singular point 

of the ^-function of the second kind, here £ = 0 is a regular 
point for all S-f unctions. Thus, for determination of the eigenvalue 
problem for £^_ = 0 in this case, also a boundary condition must be 

given. The area £p = 0 is a circular disc. If such a circular 
disc actually exists as a physical object, for instance, a circular 
screen for problems of diffraction or a circular membrane, the 
boundary condition on the disc results from the physical problem 
taken as a basis. If, however, this circular disc has geometrical 
significance only as singular surface of the coordinate system 
taken as a basis, for instance, for the determination of the acoustic 
or electrical natural oscillations inside an oblate ellipsoid of 
revolution, the eigenfunction together with its derivative must be 
required to be continuous at this circular disc which leads 
to B * 0. 



9*3 Normalization of Lame’s Wave Functions for Outside Space Problems 

One limits oneself at first to Lame's wave functions of the prolate ellipsoid of 
revolution. In the normalization of the eigenfunctions (9.2) one cannot normalize in each 
coordinate separately although the eigenfunctions contain only factors which every time are 
dependent on one coordinate only; in the element of volume 


dT = c 3 (| 2 - n 2 ) d! dn cl£p 


( 9 . 1 0 


namely, the coordinates are not separated. The condition of normalization for outside 
space problems with continuous Bpectrum reads 


nk+Ak _ f> I 
lim I dk 
S — 

Uk-Ak U 


n 1 


p2rt 


as 


an 


v 

dcp c 3 (! 2 - n 2 )u“(|,n,cp;S)*u“(S,n,9;k) 


cs '* 


Si u 


-1 u 




(9-5) 


The asterisk (*) signifies the formation of the conjugate -complex expression. 1^(7) is 

called the factor of normalization. If the indices n and m of the two eigenfunctions in. 
the integrand (9*5) would not both agree, the integral would equal zero; the same would be 
valid if the interval of integration for k would not contain the point k. However, a 
delay by proof of these properties of orthogonality is unnecessary; that proof is elementary. 

From the wave equation in elliptic coordinates one obtains the identity 



dn 


(1 - n ) 




2 Su 


dn 


u 


dn / 


(9.6) 
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For abbreviation the wave f -unction concerning the parameter values k,n,m is designated by u. 
The (^-integration of equation (9*6) gives. 2rt. The tj- integration of the second term of the 
sum on the right side of equation (9*6) gives always zero since at the ends of the interval 
of integration 1 - r\ d = 0. After carrying out the q- integration the | -integration of 
equation (9 -6) results finally in 


32 ( 1 ?- 


n 

pi 

'i2 Jt 

">1 

v 

« 

( e 2 - t ) 2 ) *1 

u*u dcp = 2 jc 

dq 

- 4* it - » ¥) 


-1 L 

0 t 

-1 



( 9 * 7 ) 


■Hie boundary condition is assumed u = finite for | = 1 in case g 1 = 1, 

or au + p — ^ = 0 in case > 1, where a and 3 shall both be real and independent 

of k. Then the content of the brackets on the right side of equation (9-7) at the 
point vanishes} for large | the asymptotic series of the 2- functions may be 

substituted. Note ihe fact that, for real boundary conditions, A/3 is always real. 

After division by S 2 - k 2 , performance of the k-integration and transition to £ = » 
there results the integral in equation (9*5) } thus one obtains as the factor of normalization 



CO 

(n + m) I V” 2 
(n - m) i /_ 2n + 2r + 1 

r=-co 


m 


m 


®n,r^)T>n,r<7) 


( 9 - 8 ) 


-a 
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ro 


The normalized wave functions are, therefore, obtained by multiplying u m (g,T},cp;k) 
by H®(y) . 

Several changes occur for the coordinates of the oblate ellipsoid of revolution. The 
eigenfunctions (9»3) have to be substituted in equation (9*5)j instead of | 2 - q 2 in 
equations (9.4), (9-5), and (9*6) one has to write g 2 + q 2 ; I 2 - 1 in equations (9.6) 
and (9*7) is to be replaced by g 2 + 1. Finally one obtains as the factor of 

normalization N^(iy) instead of N^(y) in equation (9*8) . 


9-4 Development of Lame's Wave Functions in Terms of Spherical and Cylindrical Functions 

By the following consideration one obtains a remarkable development which includes a 
number of the developments of X- and 2-functions given so far as special cases. Any Lame 
wave function can be developed in terms of such wave functions as originate by separation, 
for instance, in polar coordinates. The eigenfunction of the continuous spectrum, in 
particular, which results from equation (9-2) for B = 0 can be developed in terms of the 

/. \_m/z\ imcp 

eigenfunctions in polar coordinates \!r w ^kr;±^-je ' 0ne obtains the development 

coefficients by making, for instance, x and y — 3> 0, that is, £ — >1 and comparing 

the thus originating development with equation (5 .1) . Thus there results, if one 
expresses, moreover, x,y,z by g,q,cp and equates the coefficients with e iM P, 


n n 


00 



pss- 00 



(9-9) 

(9-10) 
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■where r in the argument. 1b given by r 2 « x 2 + y 2 + z 2 and must 
not be confused with the index r which runs .through all even 
numbers. The development (9*9) is given already by Moree (reference 6). 
One can interpret equation (9«10) aa a development of the -2- functions 
which contains still an arbitrary parameter tj. For q — >1 there 

originates, if one divides before by (l - q 2 ) m / 2 , the series (5*l)j 
for ri — >0 one obtains the series (5*42) for the W-f unctions of the 
first kind. If one differentiates equation (9*10) with respect to q 
and sets then q = 0, there results the series (5»4l) of the 
V-f unctions of the first kind. For 5 — > 80 there originates from 
equation (9*10) the series (4.1) of the X-functions of the first kind. 

At this point one can recognize why the f emulations (4.1), (4.2), 
(5*1), (5* 2 ), (5*41), and (5*42). that is, the series developments 
considered by Niven (reference 1) all had to lead to the same 
development coefficients a^.. 

Whereas equation (9*9) represents a development in terns of 
eigenfunctions in polar coordinates which have their origin at the 
point x = y = z = 0, Lame’s wave functions can be developed also 
in termB of eigenfunctions in polar coordinates with the 
origin x = y = 0, z = c . This development reads, as shown by a 
simple calculation, 


= X t (rVf(s - ’»]Ct( 6 ^i) (9.U) 

t=~00 

If one multiplies by (l - q 2 ) m ^ 2 and then sets t) = l, equation (9*11) 
is transformed into the development (5*46). For- g — > » one obtains 
a development in terns of spherical functions multiplied by sin (yq) 
and cos (yq), respectively; the special case of this development m <= 0 
is already given by Hanson (reference 5) . 

If one finally develops Lame's wave functions in terms of the 
eigenfunctions originating by separation of the wave equation in 
cylindric coordinates, there results, with the aid of equation (7.24), 


•p- 


(7) 


This integral equation vhich was derived in another way also by Kotani (reference 8) cont ains 
equation (7-24) as special case for q — >1. 

It is obvious that the developments (9*10), (9-11), and (9-12) are capable of generalization; 
one can consider complex V,|i ,7 and one can replace one or both functions of the first kind on the 
left and right side of equations (9-10), (9.11), or (9*12) by functions of the second kind; however, 
reproduction of the thus originating formulas and establishment of their domain of validity -trill be 
omitted. 


10. TEE METHOD OF GREEN'S FUNCTION FOR THE SOLUTION OF BOUNDARY 
VALUE PROBLEMS, PARTICULARLY OF RADIATION PROBLEMS 
10.1 Green's Function of the Wave Equation in Radiation Problems 


01 


io&Tj^d)^;^ \l(£ 2 _ 7 2 y :1 . n 2 )(l , a 2) 


da (9.12) 





i r T>“ 


For development of the plane wave or more generally of the spherical wave in terms of Lame 's 
wave functions of the ellipsoid of revolution one uses the method of Green's f unction thought up 
by Sommerfeld (reference 15 ) . According to this method the spherical wave is a solution of the 
inhomogeneous wave equation 


ApU(P,Q) + j£u(P,Q) = 8(P,Q) 


( 10 . 1 ) 
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the right side of which is a (recently so-called) Dirac’s 5 function. 

It was introduced first by Sommerf eld (reference 15) and designated 
by him as prong function. It has a singularity at the source 
point Q in such a manner that 

J & (P,Q) dr p = 1 (10.2) 

for each domain G which contains the source point Q, whereas 
the integral has the value zero If the domain G does not contain 
the source point. One can interpret 5(P,Q) as limiting case 
of a function which has for points of influence P in the neighbor- 
hood of the source point Q a very steep prong whereas it decreases 
toward the outside very rapidly to zero. 

The solution of equation (10.1) is for outside space problems 
uniquely determined only when besides the boundary conditions on 
the bounding areas which are at a finite distance an additional 
boundary condition at infinity is required, namely, the outgoing 
radiation condition (or else the incoming radiation condition) 
introduced by Sommerfeld (reference 15) • According to this 
condition, u(P,Q) for points P at very large distance from the 
source point Q should behave like an outgoing (or Incoming) wave. 

One designates this solution because of its special properties as 
Green's function G(P,Q;k) of the wave equation pertaining to the 
outgoing (or incoming) radiation condition. For physical reasons 
the case of the incoming radiation condition will not be considered 
below. 

All developments of this section are performed for the coordinates 
of the prolate ellipsoid of revolution^ one obtains the corresponding 
formulas for the coordinates of the oblate ellipsoid of revolution 
by replacing ^ everywhere by i 7 . 

Green’s function can be developed in the following way in terms 
of the eigenfunctions of the continuous spectrum 


G(P,Q]k) = 


drc 

k 2 - k 2 


u n(lp# 1 lpj c Ppi K ) u n^Qi r )Qi <: PQit«^ ^n(7)^ (10«3) 


n,m 


The integration over k goes from 0 to oo, the path of integration 
deviating at the point k = k in the case of the outgoing radiation 




condition into the negative- imaginary half plane. The integration over k can he performed 
according to Sommerfeld and yields for 

G(P,Qjk) = -~y~ (a“ - l^^ (3) (Cpir) e lm ( < Pp“ c PQ) ^(7) 2 [a^S^ 1 )(^7) 
n,m A L - n 

♦ b^< 2 > « Qi r)fe (1) ( V ,r) do A) 


-5 

ON 


In order to characterize the dependence of the constants A and B on n and m, the 
indices n and m were appended. For £p < £q the arguments £p and in 

equation (10.4) are to be exchanged since Green's function is symmetric in P anrl Q. 
One can recognize that there actually result outgoing waves by substituting for the 
^-functions of the third kind in equation (10.4) their asymptotic representation, since 
an outgoing wave is given by the behavior e 1 . for large r if the time dependency 
as is customary in wave physics is fixed by e" 1CDlJ . 


10.2 Development of the Spherical Wave and of the Plane 
Wave in Terms of Lame s Wave Functions 

If the space does not have any boundaries within a finite region, B » 0 in the 
eigenfunctions for the coordinates of the prolate as well as of the oblate ellipsoid of 
revolution. A may then be set equal to 1. On the other hand, it is known that Green's 
function for the entire space without boundaries within a finite region is given by the 
spherical wave 
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G(P,Q;k) 


ikr 


PQ 


-4jtr 


PQ 


(10-5) 


Therefrom one obtains immediately the development of the spherical wave in terns of Lame's wave 
functions for 


pi ij£r PQ 


?Q 


,m(3) 


Ml) 


Mi) , 




3 j£(y) (10.6) 


n.~a 


and correspondingly for < £p • Morse (reference 6) discovered this development in terms of 

another method. The development of the plane wave in terms of Lame'B wave functions originates 
from this if one moves the source point Q to infinity. For » 1 


G(P,Qjk) % 



(10.7) 


is valid where r^ is the distance of the source point from the origin of the coordinates, 

r p the radius vector of the influence point P and £q a unit vector in the direction 

of Q toward the origin of the coordinates, therefore in the direction of the direction 
of propagation of the plane wave. For equation (10.6) this limiting process is performed 

r, is used in place of German script V* -J 

e is used in place of German script 
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"by introducing for the -g-f unction of the third kind the asymptotic series (5.12) and setting 
according to equations (2.llg) and (2.11a), respectively, £q w kr^. Then there results for 

the plane .wave running in direction 


e ik(e Q ,£ p ) _ 2£ 

k 2 


^ 1 - n ^( 1 )(g p . r ) 3 ^( ;L )( t ^. y )^( 1 )(r, Q . 7 ) e im(cp Q ^ lT(r) 2 


i r b m 


n,r 


( 7 ) ( 10 . 8 ) 


n.m 


r=-oo 


The direction in which the source point is situated can he characterized instead of by the 

vector -e by the coordinates T} and cp . According to eauations (2.1g) and (2.1a) 

H Q ^ 


-(%£p) - V'-i ' 1 - -'1q 2 cos cp Q + y p \|l - ti q 2 Sin cp Q + ZpHq 

= C\f|p 2 + 1 \/( 1 - Tlp 2 )(l - TIq 2 ) COS (9q - C[>p) + c| pTJpT] q 


(10.9) 


is valid. 


10.3 Diffraction of a Scalar Spherical Wave or Plane Wave 
on the Ellipsoid of Revolution 

The method of Green's functions can also be used for treatment of problems of diffraction. 

rhi 

In the scope of a scalar diffraction theory a homogeneous boundary condition au + (3 = 0 is 

on 

prescribed on the diffracting body; at first, the problem of the determination of the eigenfunctions 
is solved. With these eigenfunctions. Green's functions (10-3) and (10.4), respectively, are 
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formed and the contention made that it solves the problem of dif- 
fraction. Actually it represents a wave which comes as spherical 
wave from the source point Q, satisfies on the diffracting surface 
the given boundary condition (as does each single term of the sum), 
and which behaves at infinity like an . outgoing spherical wave . If 
the source point Q, in particular, lies at infinity, Green's 
functions represent the superposition of a plane wave and of an 
outgoing spherical wave originating from the diffracting body with 
an amplitude which, in general, is dependent on direction. Treat- 
ment of the diffraction problem for a source point within a finite 
region is omitted. One starts immediately from equation (10.8) and 
contends that the solution of the diffraction problem of a plane wave 
at the ellipsoid of revolution is given by 


2*2 


z 


n,m 


- l £ — -*g (3) (fe»r> 


n 


m( 3) n . » 


n n ^ 


nF(y) 
n 


2 ^-n 


i r b m (y) (10.10) 


r=-oo 


in the case of the boundary condition u = 0 for £ = Under N^(y), 

one understands therein the factor of normalization (9*8) with 
M* + BB* =1. In the case of the boundary condition —I = 0 

for £ = one has to replace the two 15-functions with the 
argument £p in equation (10 .10) by their derivatives with respect 
to £p at the point The first term of the sum In the brackets 

of equation (10.10) yields, when the sum over n,m is formed, exactly 
the plane wave (10.8); the second term of the sum gives outgoing 
spherical waves; furthermore, the wave eqtxation and the boundary 
(surface) condition are satisfied by each separate term of the sum; 
the contention is therefore proved. 

For the diffraction at the infinitely thin wire of finite length, 
one has to set £ p = 1. -2,^(3) (£^.7) then becomes infinitely large 

and, In equation (10.10), there remains only the plane wave. Thus 
an infinitely thin wire does not present an obstacle for a plane wave. 
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For the diffraction at the infinitely thin circular disk, ,7 is to 
"be replaced in the formulas "by i 7 end is to he set equal to zero. 

(°;il) has a finite value so that the outgoing spherical waves 

do not disappear; that is, even an infinitely thin disk represents 
an essential disturbance for a plane wave striking it. 


11. TABLES 

11.1 Comments to the Tables 


The tables in section 11.2 contain power series developments 
to 7 10 , inclusive, for the eigenvalue 1 ^( 7 ) according to 

equation ( 6 . 13 ) and to 7 ^, inclusive, for the coefficients ( 7 ) 

and b m ( 7 ) according to equations (6 .14) to ( 6 . 21 ) and 
equation (^.11). Furthermore , to 7 °, inclusive, the coef- 
ficients a™ T f\ Q ( 7 ), according to equation ( 7 * 6 ), are given 
for all those cases where a_ 2 /a 0 , a_^/aQ, and a_g/a 0 disappear. 

As far as the values of the coefficients a r /a Q and b r /b 0 are 

not .given in the tables, they disappear; than one must use for the 
X-f unctions of the second kind the series (7*7) and the table for 
the c^/Sq. 

The region of the n- and m-values in the tables extends 
from m « 0 , 1 , 2, . . .,9 and from n = m, m + 1, . . ., 9* For 
negative m, which are integers, reference is made to the 
relation ( 4 . 12 ). 

The last given digit is, in general, probably certain; only where 
the following digit after rounding up or off, respectively, is a 5 , 
the last given digit would have to be changed in a few cases by unity. 
In the cases of the end digits ... 5 , *«50, *500/ and so forth, it 
is mostly indicated by a line over or under, respectively, the last 
digit whether the respective decimal fraction had. been originated by 
rounding up or off. 

The given broken off series developments in terms of powers 

p p 

of 7 are the more useful, the smaller y and the larger n. 
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p 

For n » 0 the series begin to be useless only at y = 10; for 
larger n they can be used up to far^ larger values of 7 2 . Below, 
a few of the first eigenvalues for y 2 = 10 are given as they 
follow from the exact numerical calculation and from the power- 
series development to 7IO, inclusive. 


n 

0 

2 

4 

^°(V io) 

2.305040 

11.790395 

25.251313 

X° (fiE) appro xima t i on 

2.215 

11.880 

25 .25147 

L_ 


Figure 1 gives a survey on the dependence of the lowest 
eigenvalues on y . 

The tables in section 11-3 are taken from the thesis of 
Bouwkcmp (reference 10) . They contain the eigenvalues X°(y) for 

a number of pairs of values n. y 2 , and ’the coefficients a m (7) 

n,r 

of the pertinent X-f unctions. These latter are fixed so that 



_ 2n + 1 
2n + 2r + 1 



= 1 


( 11 . 1 ) 


The integral of normalization then (compare equation (7. 19)) has 
the value 



( 11 . 2 ) 


These tables contain further the values X^ (lj7) and X^(0;y) 

for even and dX' ' (0 }y) / dp for odd n. The signs of the a are 

different from those of Bouwkamp since the present series (4.1) 
and (4.2) contain in the coefficients a factor i r which is missing 
in reference 10 by Bouwkcmp. 
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Since the y assume in these tables only negative values, 
these functions are appropriate for the treatment of problems 
concerning the oblate ellipsoid of revolution or for the investi- 
gation of the eigenvalues of the ion of the hydrogen molecule, 
whereas, the tables in Section 11.2 where y 2 can be positive as 
well as negative, may be used for problems of the oblate as well 
as of the prolate ellipsoid of revolution. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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11.2 — Eigenvalues 1^(7) and Development Coefficients ^^( 7 ), 


Bepresented By Broken-Off Power Series in 7 



TABLE 1 .- Jl^r) 

- 

n(n + 1 ) | X 

10 10 AS POWER 

SERIES IN 7 



m = 0, n = 0 

+ 33333353337 2 

-1181181187^ 

.+I7O3II67 6 

+ 1358687 8 

-21280. 967 10 

n = 1 

+6000000000 

— 

68571129 

- 609521 

+ 25896 


872.80 | 

n = 2 

+5238095238 

+101500918 

-1760812 

-IIIO89 

+21112.17 

n = 3 

+5111111111 


32911763 

+ 595989 

- 26512 

— 

887.76 

n = U 

+5061-935062 

+ 

17750507 

+ 53117 

+ 5010 

- 

132.19 

n * 5 

+ 501273501-3 

+ 

II298966 

+ 11653 

• * + 577 

+ 

11.80 

n = 6 

+5030303030 

+ 

787113! 

+ 365 S 

+ 150 

+ 

0.91 

n = 7 

+502262II 31 

+ 

5819121 

+ 1113 

+ 53 

+ 

0.1s 

n = 8 

+5017513860 


II8265I - 

+ 628 

+ 22 

+ 

0.03 

n = 9 

+5011005602 

+ 

3562110 

+ 310 

+ 11 

+ 

0.01 

m = 1, n = 1 

+2000000000r 2 

_ 

I571I2867 4 

+ 12190187 6 

- 2IO3I7 8 


205. 717 10 

n = 2 

+1285711286 

— 

38872692 

+ 111210 

+ 5682 

— 

76.28 

n = 3 

+466666666jT 

+ 

13617587 

-II825OI 

+ 21357 

+ 

229.72 

n = h 

+1805191805 

+ 

11902117 

- 131503 

- 5669 

+ 

78.30 

n a 5 

+1871791872 

+ 

889I60I 

- 31166 

- 337 

— 

23.6s 

n = 6 

+1909090909 

+ 

6698782 

- 10150 

26 

_ 

1.93 

n = 7 

+1932126697 

+ 

517^730 

- 1008 

+ 6 

_ 

0.31 

n = 8 

+1917368121 

+ 

I099I03 

- 1807 

+ 7 

— 

0.08 

n = 9 

+1957983193 

+ 

3320032 . 

899 - 

+ 5 

- 

0.02 

1 = 2, n = 2 

+H285711297 2 

_ 

19I363I67 1 * 

+ 3606007 6 

- 58227 s 

+ 

57 . 717 10 

n = 3 

+3333333333 

- 

2211668 9 

+ 127901 

+ 607 

— 

22.60 

n = U 

+1025971026 

— 

213987I 

- 309653 

+ 6053 

— 

56.26 

n = 5 

+1358971352 

•f 

258I89I 

- 101857 

- 532 

+ 

23.27 

n = 6 

+1515151515 

+ 

3I76672 

- 39IOI 

- 205 

— 

1.22 

n = 7 

+1660633181 

+ 

3361211 

- 16771 

67 

— 

O.58 

n = 8 

+1736812105 

+ 

3005596 

- 7911 

22 

— 

0.30 

n = 9 

+1789915966 

+ 

262O83I 

- 1052 

8 

- 

0.07 

m = 3, n = 3 

+ 11111111117 2 

— 

99763067^ 

+ 1326387 s 

- 16987 8 

+ 

17 . 807 10 

n = U 

+2727272727 

— 

13870127 

+ 76121 

83 

— 

1.51 

n = 5 

+3501273501 

— 

I9200I0 

- 92311 

+ 1750 

— 

18.22 

n = 6 

+3939393939 

— 

877352 

- 51168 

+ 119 

+ 

1.55 

n = 7 

+I208HI796 

+ 

755560 

- 27781 

32 

+ 

0.53 

n = 8 

+I38596I912 

+ 

1369050 

- ill 31 

27 

— 

0.01 

n = 9 

+I509803922 

+ 

15I8920 

- 7813 

11 


0.01 

m = 4 , n = 4 

+ 9090909097 2 

_ 

5779 

+ 573167 6 

- 5787 ® 

+ 

5 . 217 10 

n = 5 

+2307692308 

- 

9I03323 

+ II360 

- 118 

— 

0.79 

n = 6 

+3090909091 

— 

1839057 

- 28865 

+ 570 

— 

5.53 

n = 7 

+3571660633 

— 

2037903 

- 26567 

+ 129 

+ 

0.73 

n - 8 

+389I7368I2 

— 

5305^0 

- 17172 

+ 17 

+ 

0.29 

n = 9 

+U1761 70 59 

+ 

211118 

- IOI67 

1 

+ 

0.06 

m = 5 , n = 5 

+ 7692307697 2 

_ 

36H3297 4 

+ 278837 s 

— 2II7 8 

+ 

1 . 707 10 

n = 6 

+2000000000 

— 

627I510 

+ 26H9 

80 

— 

0.07 

n = 7 

+2760180995 

- 

1157529 

- 8311 

+ 205 

- 

1.82 

n = 8 

+3263157892 

- 

2301599 

- 1278s 

+ 79 

+ 

0.03 

n = 9 ! 

+3613II5378 

- 

1096I9S 

- 10175 

+ 22 

+ 

0.10 

m = 6, n = 6 

+ 6666666677 ^ 



2HOO877 4 

+ 1I8IO7 6 

977 ® 

+ 

0.627 10 

n = 7 

+1761 705882 

- 

1199 311 

+ 16310 

19 

+ 

0.01 

n = 8 

+2191228070 

— 

3II0673 

- 1336 

+ 80 

— 

0.65 

n = 9 

+2997198880 


2221592 

- 6051 

+ 11 

- 

0.07 

m = 7 , n = 7 

+ 5882352967 2 


17110357 4 

+ 8 l 737 S 

167 ® 

+ 

0 . 2 S 7 10 

n = 8 

+1578917368 

— 

3332131 

+ 10135 

29 

+ 

0 .C 5 

n = 9 

+2268907563 

- 

2822595 

+ 967 

+ 32 

- 

0.25 

m = 8, n = 8 

+ 5263157897 2 

— 

12I96627 4 

+ 51177 S 

237 ® 

+ 

0 . 117 10 

n = 9 

+1I28571I-29 

- 

2535176 

+ 6898 

18 

+ 

0.03 

m = 9 , n = 9 

+ I76 190I 767 s 

- 

9 3895^7 4 

+ 32367 s 

137® 

+ 

0 . 027 10 
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TABLE 2.- 


(7) 

w 


x io 10 


AS POWER SERIES IN 


7 


m.= 0, n = 0 

+ 11111111117 2 

- 35273 3697 1 * 

- 10100087 6 

n = 1 

+ 400000000 

+ 3555556 

— 

151062 

n = 2 

+ 244897959 

+ 30290)1 

+ 

208818 

n = 3 

+ 17636681+3 

+ 66996 

+ 

1*6389 

n = b 

+ 13774101+7 

+ 21683 

+ 

18710 

n = 5 

+ II2963959 

+ 8738 

+ 

91+67 

n = 6 

+ 9 57261+96 

+ 1+071 

+ 

5462 

n = 7 

+ 8301+1+983 

+ 2105 

+ 

3439 

1 n = 8 

+ 73325729 

+ H 77 

+ 

2306 

n = 9 

+ 6561+0291 

+ 700 

+ 

1622 

m = 1, n = 1 

+ 1333333337 2 

- J 5555567 U 

+ 

6131+87^ 

n = 2 

+ 1221+1+8980 

- U 51+356 

— 

17898 

n > 3 

+ 105820106 

- 120593 

+ 

16210 

n = b 

+ 91827362 

- 1+3366 

+ 

9870 

n - 5 

+ 8068851+2 

- 18723 

+ 

5930 

n = 6 

+ 71791+872 

9160 

+ 

3767 

n = 7 

+ 61+59051+2 

1+912 

+ 

2524 

n = 8 

+ 58660583 

2826 

+ 

1768 

n = 9 

+ 53705693 

1719 

+ 

1285 

m = 2, n = 2 

+ I08163277 2 

- 7572607 ^ 

+ 

122277 6 

n =+ 3 

+ 52910053 

- 301482 

— 

1430 

n = 1+ 

+ 550961+19 

- 130098 

+ 

2386 

n = 5 

+ 53792361 

- 621+11 

+ 

2541 

n = 6 

+ 51282051 

- 32716 

+ 

2057 

n = 7 

+ 1+81+1+2907 

- 181+20 

+ 

1579 

n = 8 

+ 1+5621+898 

- 10990 

+ 

1207 

n = 9 

+ 1+2961+551+ 

6877 

+ 

931 

ra • 3 , n = 3 

+ I763668H7 2 

- 23!+ 1+867** 

+ 

30037 6 

n = 4 

+ 2751*8209 

- 151781 

+ 

165 

n = 5 

+ 322751+17 

- 87376 

+ 

621 

n = 6 

+ 31+188031+ 

- 50891 

.+ 

828 

n = 7 

+ 31*602076 

30700 

+ 

812 

n = 8 

+ 31+218671+ 

- 19232 

+ 

718 

n = 9 

+ 331+16876 

- 121+81 

+ 

609 

m = U, n = U 

+ 91827367 2 

- 910687*1 

+ 

9187 6 

n = 5 

+ 16137708 

- 78638 

+ 

210 

n = 6 

+ 20512821 

- 51+963 

+ 

242 

n = 7 

+ 23068051 

3681+0 

+ 

329 

n = 8 

+ 21+1+1+1910 

- 21+727 

+ 

361 

n = 9 

+ 25062657 

— 168^9 

+ 

353 

m = 5 , n = 5 

+ 53 T 92367 2 

- 411927 ** 

+ 

3327 6 

• n = 6 

+ 102561+10 

- 1+3182 

+ 

131 

n = 7 

+ 138I+O83O 

- 34735 

+ 

118 

n = 8 

+ 16291+606 

- 25905 

+ 

152 

n = 9 

+ 17901898 

- 18912 

+ 

178 

m = 6, n = 6 

+ 3I+I88037 2 

- 207937** 

+ 

1367 6 

n = 7 

+ 69201+15 

- 25086 


75 

n = 8 

+ 9776761 * 

221+51 

+ 

64 

n = 9 

+ 11931+598 

18212 

+ 

79 

m = 7 , n = 7 

+ 27068057 2 

111 + 037 ** 

+ 

627 6 

n = 8 

+ 1+888382 

15307 

+ 

1 : 

n = 9 

+ 7160TO 

- 11+900 

+ 

37 

i = 8, n = 8 

+ 1629U6I7 2 

66727*1 

+ 

307 6 

n = 9 

+ 3580380 

9 - 71+3 

+ 

25 

m = 9 , n = 9 

+ LL93I16O7 2 

1 + 111 + 7 ** 

+ 

167 6 




NACA TM No. 122^ 


a n li^) in 
TABIZ 3.- -£*- x 10 10 


%6 (r) i 

AMD X 10 1 


as povxs sir iks 




“4/ a o 

x 10 10 


«6/>o x 10 ' 

- 0 

+190476197* 

■ -7.696OI7 6 

+I374297 6 

= 1 

+ 4535147 

+ 

55817 

+ 24667 

• 2 

+ 2061431 

+ 

3740 

+ 8970 

- 3 

+ II77274 

+ 

684 

+ 4262 

- 4 

+ 760702 

+ 

189 

+ 2355 

- 5 

+ 53159S 

+ 

67 

+ 1437 

- 6 

+ 392250 

+ 

28 

+ 940 

- 7 

+ 301245 

+ 238563 

+ 

13 

+ 649 

* 8 

+ 

7 

+ 466 

- 9 

+ 193571 

+ 

4 

+ 346 


+ 9070297 4 
+ 687144 
+ 504546 
+ 380350 
+ 295331 
+ 235350 
+ 191701 

+ 159042 

+ 134010 

+ 1374897 * 
+ 168182 
+ 163002 
+ 147663 
+ 130750 
+ 115021 
+ 101209 


334907 

3740 

879 

284 

111 

52 

24 

13 

7 


+ 35247° 
+ 2242 
+ 1421 
+ 942 


89340 

- 

2 S 

+ 

119 

336367* 

_ 

684 7 6 

+ 

517 

54336 

- 

473 

+ 

79 

6328s 

— 

277 

+ 

87 

65375 

— 

l6l 

+ 

85 

63900 

— 

95 

+ 

79 

60725 

— 

58 

+ 

72 

56853 


37 

+ 

64 

108627* 

_ 

17O7 6 

+ 

117 

21095 

— 

166 

+ 

22 

28013 

- 

124 

+ 

28 

31950 

— 

86 

+ 

32 

33736 

— 

58 

+ 

33 

34112 

- 

40 

+ 

32 

42197* 

_ 

527 6 

+ 

37 

9339 

— 

6S 

+ 

7 

13693 

— 

58 

+ 

11 

16868 

— 

46 

+ 

13 

18951 

“ 

3S 


IS 

I8687* 

_ 

197 6 

+ 

17 

4564 

— 

28 

+ 

3 

7222 

— 

28 

+ 

4 

9475 

- 

23 

+ 

6 


m = 9, n = 9 


+ 


+ 

+ 
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TABLE 4.- 


3r-g (y) x IQ 10 


,(r) 


AS FOWXR SKRIE3 IX 


r 


a * 0, a = 2 

- 2222222227 2 

+ 

70546747^ 

+ 

1519717 6 

n = 3 

- 171428571 

— 

1523810 

+ 

42927 

n = 4 

- 136054422 

— 

168280 

— 

26152 

a = 5 

- 112233446 

— 

42634 

— 

8241 

a = 6 

- 95359186 

- 

15011 

- 

4452 

n = 7 

- 82840237 

- 

6408 

- 

2790 

n a 8 

- 73202614 

- 

3113 

— 

1889 

n = 9 

- 65561828 

_ 

1662 

- 

1346 

a = 1, a = 1 

-3333333333r 2 

-13333333337 1 * 

-540952 38I7 6 

a » 2 

— 666666667 

- 

63492063 

- 

6478782 

a = 3 

- 3^2857143 

+ 

9142857 

- 

197018 

n * k 

- 226757370 

+ 

841400 

+ 

16156 

n - 5 

- 168350168 

+ 

191852 

- 

7170 

n * 6 

- 133502861 

+ 

63047 

— 

5138 

n - 7 

- 110453649 


25630 

- 

33«I 

n * 8 

- 94117647 

+ 

12009 

— 

2295 

n . 9 

- 81952286 

+ 

6232 

- 

1621 

a = 2. n = 2 

-13333333337 2 

— 

6349206357^ 

_ 

8055289O7 6 

a » 3 

- 571428571 

+ 

76190476 

+ 

8760622 

n = k 

- 340136054 

+ 

6310502 

— 

122277 

a =* 5 

- 235690236 

+ 

1342964 

— 

4268 

n * 6 

- 178003814 

+ 

420316 

- 

4295 

a - 7 

- 142011834 

+ 

164766 

— 

3369 

a = 8 

- 117647059 

+ 

75054 

- 

2449 

a * 9 

- 100163905 


38087 

- 

1778 

a » 3» a » 3 

- 8571428577 2 

+ 

2666666677k 

+ 

8838009^7^ 

n * 4 

- 476190476 

+ 

20614306 

— 

2451288 

n - 5 

- 3H*25 3648 

+ 

4178111 

— 

64139 

n * 6 

- 228862047 

+ 

1260948 

- 

7860 

a - 7 

- 177514793 

+ 

480568 

— 

3572 

n = 8 

- 14379085a 

+ 

214043 

— 

2475 

n = 9 

- 120196685 

+ 

106642 


1829 

a = 4, a = 4 

- 6 349206 357 2 

+ 

494743357k 

_ 

129516487 6 

a = 5 

_ 404040404 

+ 

9669343 

— 

380104 

n = 6 

- 286077559 

+ 

2837133 

— 

34429 

a = T 

- 216962523 

+ 

1057249 

- 

6853 

a = 8 

- 172549020 

+ 

462332 

— 

2950 

a > 9 

- 142050628 

+ 

226858 


1908 

» » 5, a - 5 

- 5050505057 2 

+ 

189933527k 

_ 

136583I7 6 

a « 6 

- 349650350 

+ 

5449096 

— 

123724 

a » T 

- 260355030 

+ 

1993669 

— 

19430 

n » 8 

- 203921569 

+ 

858612 

- 

5180 

a » 9 

- 165725733 

+ 

415905 

— 

2347 

a = 6 . a - 6 

- 4195804207 s 

+ 

94451007k 

_ 

347OH7 6 

a - 7 

- 307692308 

♦ 

3403334 

— 

52456 

a - 8 

- 237908497 


1446929 

— 

11529 

a » 9 

- 191222000 

+ 

693176 


3762 

a =* 7, a =. 7 

- 3589743597 2 

+ 

54143957k 

_ 

1238787 6 

n * 8 

- 274509804 

+ 

2276636 

— 

25809 

a - 9 

- 218539428 

+ 

1080274 

“ 

7154 

a * 8. a * 8 

- 313T254907 2 


34024467k 


537207 6 

a - 9 

- 247678019 

♦ 

1601021 

- 

14022 

a - 9, a - 9 

- 2786377717 2 

+ 

22814567k 

- 

265037 6 
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TABU 5.- 


S,o ( ^ 


x 10 


10 


AKD 


\o<^ 


X 10 10 


AS FCWXR SEtHS DT 


7 



a -4/*o 

X 10 10 


“*-6/ a o x 1q1 ° 

a * 0, n = 4 

+ 9070297^ 


67317 6 


n = 5 

+ 68711+1* 

+ 

2610 


n * 6 

+ 5045lt6 

+ 

371 

16027 6 

n = 7 

+ 380350 

+ 

106 

1207 

n * 8 

+ 295331 

+ 

39 

852 

n = 9 

+ 235350 

+ 

17 

610 

* * lj a * 3 

+ 1901*76197* 

+ 

33B6247 6 


n = U- 

+ 453511*7 

+ 

100968 


n > 5 

+ 20611*31 

— 

23492 

- 458107° 

n * 6 

+ 1177271* 

— 

2595 

- 11212 

n - 7 

+ 760702 

— 

635 

4830 

n = 8 

+ 531595. 

— 

213 

2557 

n « 9 

+ 392250 

- 

86 

1524 

i ■ a, n . 2 

-666666667-** 

-1269841277 6 


n = 3 

+ 95238095 

+ 

8465608 


n = 

+ 136051*1*2 

+ 

1514520 

+ 15117167° 

n = 5 

+ 4810005 

- 

274074 

- 320667 

n * 6 

+ 2354548 

— 

25945 

- 44849 

n - 7 

+ 1369260 

— 

5719 

- 14490 

n * 8 

+ 885992 

— 

1776 

6392 

n = 9 

+ 616393 

- 

677 

3353 

* * 3. n * 3 

+2857142867* 

+ 592592597 6 

-317460327 6 

n = 4 

+ 31746032 

+ 

8245723 

+10582011 

n - 5 

+ 9620010 

— 

1279014 

- 1234668 

n = 6 

+ 4238186 

- 

108971 

- 134546 

n - 7 

+ 2282100 

— 

22241 

- 36224 

n ■ 8 

+ 1392273 

- 

6514 

- 14063 

n - 9 

+ 924590 

— 

2370 

6705 

a - 4, n <• 4 

Jl 

+ 634920637 

♦ 

296846017 6 

+423280427 6 

iv ■ ■ 5 

+ 17316017 

— 

4144004 

- 3848004 

n * 6 

+ 7063643 

— 

326912 

- 336364 

n - 7 

+ 3586157 

- 

62911 

- 79692 

n ■ 8 

+ 2088409 

- 

17587 

- 28127 

n . 9 

+ 1335519 


6162 

- 12452 

a » 5, n - 5 

+ 288600297* 


108533447 6 

- 962OOIO7 6 

n ■ 6 

+ 11100011 

— 

807274 

- 740001 

n « 7 

+ 5379236 

— 

148289 

- 159385 

n » 8 

+ 3016591 

- 

39912. 

- 52235 

n - 9 

+ 1869726 


13555 

- 21792 

a = 6 , n » 6 

+ 166500177* 

_ 

174909 37 6 

- 14800017 6 

n - 7 

+ 7770008 

— 

309394 

- 296000 

| n * 8 

+ 4223228 

— 

80721 

91412 

n « 9 

+ 2549627 

- 

26700 

- 36319 

a - 7, n = 7 

+ 108780117* 


59066I7 6 

- 5180017 6 

n * 8 

+ 5758947 

- 

150108 

- 152353 

n = 9 

+ 3399502 

” 

48546 

- 58111 

a - 8, n = 8 

+ 77102387* 


2628057 6 

- 2447697^ 

n = 9 

+ 4445503 

- 

83016 

- 89808 

m = 9, n = 9 

+ 57156477* 

- 

1351977 6 

- 1347127 6 
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TABLE 6 .- NUMERICAL VALDES OF THE COEFFICIENTS oj r /aJ Q 
APPEARING IN THE SERIES DEVELOPMENTS (j.j) 



a -kj & o X 10 10 

E : 0, n * 0 

+1666666677^ 

-246913587 6 

n = 1 

-277777778 

+31746032 

n * 2 

+ 55555556 

- 705467 

n = 3 

+ 4761905 

- 28219 

at = 1 , n = 1 

+27777777874 

+952380957 6 

n = 2 

—I 66666667 

- 6349206 
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[jor 

2 ( *> ln 
wnTf 7 v 

AS POWER SERIES IB 7 
compar® tharafore tabla 

1 1 

N 

" (r) 

n,o 

a » 0, th«r* 1 b t> r » a r j 

m » 1, n 

- 1 

+8OOOOOOOO7 2 

-213333337 4 

+3680897 6 

n 

= 2 

+408163265 

- 1514520 

- 59660 

n 

- 3 

+264550265 

- 301482 

+ 40526 

n 

- 4 

+192837U66 

- 91068 

+ 20727 

n 

- 5 

+150618612 

- 3^950 

+ 11069 

n 

= 6 

+123076923 

- 15704 

+. 6458 

n 

= 7 

+103806228 

789^ 

+ 4057 

n 

=i 8 

+ 89620336 

4357 

+ 2702 

n 

= 9 

+ 787683^9 

2522 

+ 1881 

a = 2, a 

= 2 

+61224489 87 s 

-1135890 37 4 

+1834047 6 

n 

= 3 

+370370370 

- 2110372 

- 10012 

n 

= 4 

+257116621 

- 607123 

+ 11136 

n 

= 5 

+193652501 

- 224681 

+ 9149 

n 

= 6 

+153846154 

- 98147 

+ 6171 

n 

= 7 

+126874279 

48243 

+ ^135 

n 

« 8 

+107544403 

- 25901 

+ 2844 

n 

= 9 

+ 93089868 

- 14900 

+ 2017 

m = 3, n 

= 3 

+49382716O7 2 

- 656560 37^ 

+ 840717 6 

n 

= 4 

+330578512 

- 1821369 

+ 1978 

n 

* 5 

+242065627 

- 655319 

+ 4654 

n 

= 6 

+188034188 

- 279902 

+ 4553 

n 

= 7 

+152249131 

- 135080 

+ 357*+ 

n 

= 8 

+127097931 

- 71^32 

+ 2666 

n 

= 9 

+108604845 

40562 

+ 1980 

m = 4, n 

= 4 

+4132231407 2 

- 409 80817 ^ 

+ 412937 6 

n 

=■ 5 

+295857988 

- 1441703 

+ 3851 

n 

= 6 

+225641026 

- 604588 

+ 2666 

n 

= 7 

+179930796 

- 287353 

+ 2567 

n 

= 8 

+148280919 

- 150013 

+ 2193 

n 

= 9 

+125313283 

- 84244 

+ 1766 

m = 5, n 

- 5 

+355029586r 2 

- 27186397 ^ 

+ 2188O7 6 

n 

= 6 

+266666667 

- 1122807 

+ 3410 

n 

= 7 

+209919262 

- 526814 

+ 1790 

n 

=* 8 

+171093368 

- 272002 

+ 1601 

n 

-9 

+143215181 

- 151296 

+ 1422 

n • 6, a 

- 6 

+3iiiiini7 2 

- 18921387^ 

+ 123857 6 

n 

- 7 

+242214533 

- 878023 

+ 2623 

n 

* 8 

+195535278 

- 449019 

+ 128? 

n 

- 9 

+162310538 

- 247677 

+ 1072 

» = 7, a 

- 7 

+2786166097 2 

- 13683477 4 

+ 741O7 6 

n 

=* 8 

+221606648 

- 693939 

+ 1939 

n 

- 9 

+I82599356 

- 379959 

+ 952 

a « 8, a 

= 8 

+2493074797 2 

- 10208917^ 

+ 46437 6 

n 

- 9 

+204081633 

- 555324 

+ 1422 

a - 9 , a 

- 9 

+2267573707 2 

- 7815677 ^ 

+ 3025V 6 
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table 8.- x io 10 


AID 


K 0 ^) 

n,o 


x io 10 

O r) 


AS P OWR SIR US U 7 



B 


0, there le bj. 


a r j compare therefore table 3 


3 



bjb 0 x IO 10 

b 6 /b 0 X IO 10 

b » 1, n » 1 

+ 136051127 ^ 

- 5023527 6 

+ 9866 T 7 6 

n => 2 

+ 1810002 

- 26177 

+26909 

n » 3 

+ 2351518 

- 1101 

+IO65I 

n = U 

+ 1369260 

- 1021 

+ 5181 

n - 5 

+ 885992 

- 333 

+ 2871 

n a 6 

+ 616393 

- 130 

+ 1716 

n - 7 

+ 151867 

- 58 

+ 1135 

n - 8 

+ 311592 

28 

♦ 777 

n - 9 

+ 27099s 


+ 551 

B a 2 , n a 2 

+ 96200107“ 

-2617697 6 

+ 672737 6 

n a 3 

+ 42 38186 

- 36931 

+23I39 

n *= U 

+ 2282100 

- 8508 

+10362 

n = 5 

+ 1392273 

- 2615 

+ 5337 

n a 6 

+ 92 U 590 

- 972 

+ 3056 

n - 7 

+ 652697 

- 117 

+ 1892 

n « 8 

+ 182128 

- 198 

+ 1213 

n a 9 

+ 3695UI1 

- 102 

+ 856 

m a 3, n a 3 

+ 70636137“ 

-H36337 6 

+1687876 

n a h 

+ 3586157 

- 31198 

+ 19212 

n - 5 

+ 2088^09 

- 9151 

+ 9339 

n * 6 

+ 1335519 

- 3282 

+ 5091 

n = 7 

+ 913776 

- 13 fo 

+ 3027 

n = 8 

+ 657857 

630 

+ 1921 

a a 9 

+ 1*92725 

- 317 

+ 1281 

b a I*, n = it 

+ 5 3792 367^ 

- 812 337 6 

+ 336787^ 

n a 5 

+ 3016591 

- 23800 

+15565 

n = 6 

+ 1869726 

- 8277 

+ 8150 

n = 7 

+ 1216058 

- 3313 

+ !6tS 

n = 8 

+ 8771 *t 2 

- 1512 

+ 2881 

n = 9 

+ 6UI4.333 

- 7 l 7 

+ 1876 

m a 5, n = 5 

+ 12232287 ** 

- 5235 S 7 6 

+219027 6 

n - 6 

+ 25 lt 9627 

- 17737 

+12592 

n = 7 

+ 166 lit 11 

- 7002 

+ 7017 

n = 8 

+ lib 70 32 

- 310^ 

+ 1211 

n = 9 

+ 828128 

- 1509 

+ 2680 

m a 6, n - 6 

+ 3 3995027 ** 

- 3 H 597 6 

+188927 6 

n = 7 

+ 2172611 

- 13231 

+10256 

n = 8 

+ H7I756 

- 5770 

+ 6016 

n = 9 

+ IOU93I2 

- 2761 

+ 3752 

m » 7 , n = 7 

+ 27933617“ 

- 231977 6 

+ H65I7 6 

n =* 8 

+ 1868021 

- 9966 

+ 8123 

n a 9 

+ 1311678 

- 4706 

+ 5160 

B = 8, n = 8 

+ 23350307 U 

- 1629O7 6 

+ 115817 6 

n a 9 

+ 1620308 

- 7602 

+ 6981 

n = 9 j n * 9 

+ 19803767“ 

- II7697 6 

+ 930876 
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TABLE 9.- 


x ioio 

b m 
n,o 


(7) 


AS POWER SERIES IN 


7 



m = 0, there Is 


b r = a r ; 


compare therefore table 




m = 

1 , 

n 

= 3 

-571428577 2 

+I5238IO7 14 ' 

-3283676 



n 

= 4 

-68027211 

+ 252420 

+ 4847 



n 

= 5 

-67340067 

+ 76741 

- 2868 



n 

= 6 

-63572791 

+ 30022 

- 2447 



n 

= 7 

-59171598 

+ 13730 

- 1811 



n 

= 8 

-54901961 

+ 7005 

- 1339 



n 

= 9 

-50992534 

+ 3878 

- 1009 

m = 

2, 

n 

» 4 

-226757377 2 

+ 4207007^ 

- 815276 



n 

= 5 

-33670034 

+ 191852 

- 610 



n 

= 6 

-38143674 

+ 90068 

- 920 



n 

= 7 

-39447732 

+ 45768 

- 936 



n 

= 8 

-39215686 

+ 25018 

- 816 



n 

= 9 

-38244400 

+ 14542 

- 679 

m = 

3, 

n 

= 5 

-1122334^7 2 

+ 1492187 24 

- 229176 



n 

= 6 

-19071837 

+ IO5079 

- 655 



n 

= 7 

-23668639 

+ 64076 

- 476 



n 

= 8 

-26143791 

+ 38917 

- 453 



n 

= 9 

-27317428 

+ 24237 

- 4 16 

m = 

4, 

n 

= 6 

- 63572797 2 

+ 6 30477 ^ 

- 7657 6 



n 

- 7 

-II83432O 

+ 57668 

- 374 



n 

= 8 

-15686272 

+ 42030 

- 268 



n 

= 9 

-18211619 

+ 29084 

- 24^ 

m = 

5, 

n 

= 7 

- 39447737 2 

+ 3020Jy k 

- 29476 



n 

= 8 

- 7843137 

+ 33024 

- 199 



n 

= 9 

-IO92697I 

+ 27422 

- 152 

m = 

6, 

n 

= 8 

- 26143797 2 

+ I59OO7 4 

- 12776 



n 

= 9 

— 5463486 

+ 19805 

- 107 

m = 

7, 

n 

= .9 

- I82II627 2 

+ 90027^ 

- 6076 


For n = m and n = m + 1, b_2/b Q disappears 
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TABLE 10.- 




X 10 10 AND ll 

£ A7) 


X 10 10 AS 


PCWER SERIES IN 7 


[For m = 0, there is b r = a^, ; compare therefore table 5] 



*uK x 1010 

h - 6 l b o x lol ° 

m = 1, n = 5 
n = 6 

n = T 
n = 8 
n = 9 

+1374297^ 

+168182 

+163007 

+14766*5 

+130755 

-I5667 6 

- 371 

- 136 

- 59 

- 29 

_17276 

-213 

-203 

m = 2, n = 6 
n = 7 
n = 8 
n = 9 

+ 336367^ 

+ 54336 
+ 63285 
+ 65375 

-3717 6 
-227 
-127 
- 72 

~ 307 6 
~ 51 

m = 3, n = 7 
n - 8 

n = 9 

+ 108677 i+ 
+ 21095 
+ 28018 

-IO67 6 

- 99 

- 72 

- 77 6 

m = 4, n = 8 
n = 9 

+ 421 97^ 
+ 9339 

- 36 7 6 

- 43 


m = 5, n = 9 

+ 18687^ 

— 1476 



For n = m, n = m + 1, n = m + 2, n = m+3. 


t-kfio 


disappears . 


For n=m, n=i+l, n=a+2, n=i+3, n=jn+l(, 
n = m + 5, b_g/b 0 disappears. 
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11.3 — Numerical Magnitude of the Eigenvalues find the Development 
Coefficients for Different n,7 and m = 0 

TABLE 11 .- EIGENVALUES X°(7) 

n 


-7 2 

X 0 

X 1 

X.g 

X 3 

3 

- 1.144334 

+ 0.140119 

+ 4.530790 

+IO.494513 

4 

- 1.594507 

- .505244 

+ 4.091201 

+10 .003864 

5 

- 2.079939 

- 1.162422 

+ 3.677958 

+ 9 .517981 

6 

- 2.599717 

- 1.831051 

+ 3.288927 

+ 9.036338 

7 

- 3.151917 

- 2 .510421 

+ 2.923314 

+ 8.558395 

8 

- 3.734090 

- 3.200050 

+ 2.578205 

+ 8.O83615 

9 

- 4.343439 

- 3.899400 

+ 2.250704 

+ 7-611465 

10 

- 4.976896 

- 4.607952 

+ 1.938379 

+ 7.141428 

15 

- 8.42084 

- 8.27180 

+ .49949 

+ 4.80616 

20 

-12 .16294 

-12 .09943 

- .91127 

+ 2.45867 

25 

-16 .07904 

-16 .05041 

- 2.44860 

+ .06093 

50 

-36 .90015 

-36 .89912 

-13 .56548 

-13 .21675 

100 

-81 .02794 

-81 .02794 

-45 .48967 

-45 .48391 


2 

-7 

X4 

X 5 

x 6 

X 7 

4 


+28.00092 



5 

+17.511596 


+39 .504499 


10 

+15 .11342 

+25 .06949 

+37 .04822 


15 

+12 .81726 

+22 .68771 

+34 .63123 


20 

+10.64634 

+20.36028 

+32 .25386 


25 

+ 8.63040 

+18.08457 

+29 .91689 


50 

+ .94568 

+ 7.25075 

+18 .92267 


100 

-16 .06556 

-15.32812 

+ 2.57368 

+11.45564 


- 7 2 

x 8 


X 10 

X 12 

25 

+59 .736180 


+97.652659 

+143 .606898 

100 

+26 .56408 

+43 .49374 

+62 .82728 
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TABLE 12.- NUMERICAL VALUES OF THE COEFFICIENTS a° (7), AND 

o,r 

THE FUNCTION VALUES X^(l;r) AND X^(0;r) 


r 

v 2 -3 

-7 =3 

4 

5 

6 

0 

+0 .987210 

+O.976788 

+O.963507 

+0 .947840 

2 

- .356220 

- .478334 

- .597277 

- .710545 

4 

+ 18683 

+ 33565 

+ 52482 

+ 74937 

6 

408 

979 

1914 

3279 

8 

+ 5 

+ 16 

+ 38 

+ 79 

10 




1 

XoU) 

+1.362526 

+1.489682 

+1.615218 

+1.736681 

*o(°) 

+ .815979 

+ .74990 6 

+ .683961 

+ .619666 


r 

2 _ 
-7 =7 

8 

9 

10 

0 

+0.930429 

+0 .911986 

+0.892960 

+0 .874035 

2 

- .816037 

- .912632 

- .999698 

-l .077418 

4 

+ .100273 

+ .127817 

+ .156894 

+ .186943 

6 

5114 

7440 

- 10250 

- 13535 

8 

+ 144 

+ 238 

+ 369 

+ 541 

10 

3 

5 

8 

14 

Xo(l) 

+1.852000 

+I.96OH8 

+2.060379 

+2 .152486 

x 0 <°) 

+ .558452 

+ .501339 

+ .448846 

+ .401345 


r 

-7 2 = 15 

20 

25 

50 

100 

0 

40. 78915 

+0.72576 

+0 .67909 

+0.55601 

+0.46036 

2 

-1.34978 

-1.49587 

-1.57800 

-1 .68750 

-1 .64674 

4 

+ .33881 

+ .47816 

+ .60025 

+1.01470 

+1.38412 

6 

- 3608 

- 6622 

- .10110 

- .29658 

- .63255 

8 

+ 260 

+ 515 

+ 965 

+ 5161 

+ .18483 

10 

10 

26 

59 

595 

- 3745 

12 

+ 

+ 1 

+ 3 

+ 48 

+ 555 

14 




3 

63 

16 





+ 1 


+2 .51654 

+2.77142 

+2 .96871 

+3 .61286 

+4.35229 

x Q (°) 

+ .23073 

+ .13779 

+ .08609 

+ .01284 

+ .00081 
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TABLE 13 .-r NUMERICAL VALUES OF THE COEFFICIENTS a? (7), AND 

i,r 

THE FUNCTION VALUES AND 


r 

-7 2 = 3 

1 

5 

6 

0 

+O.997105 

+0 .991981 

+O.99238O 

+0.989330 

2 

- .116098 

- .152711 

- .188018 

- .222236 

4 

+ 3902 

+ 6812 

+ 9751 

+ 11717 

6 

63 

150 

262 

I73 

8 

+ 1 

+ 2 

+ 1 

+ 9 

Xl(l) 

+I.II7169 

+1.15H659 

+I.190H5 

+1 .226765 

Xl(O)' 

+ .830138 

+ .778367. 

+ .728028 

+ .682557 


r 

V 5 = 7 

8 

9 

10 

0 

+O.9859IO 

+0 .982170 

+0.978150 

+0 .973908 

2 

- .255039 

- .286170 

— 

.316613 

- .315385 

1 

+ 19595 

+ 25008 

+ 

30907 

+ 37211 

6 

753 

- 10 66 

— 

1177 

1969 

8 

+ 16 

+ 27 

+ 

12 

+ 62 

10 



- 

1 

1 

x^d 

+1.261293 

+1.29I7I1 

+1 

.327190 

+1.358536 

Xi(0)» 

+ .638528 

+ .597088 

+ 

.558051 

+ .521111 


r 

E 

2 = 15 

20 

25 

50 

100 

0 

+0.95067 

+0 .92651 

+0.90339 

+0.81333 

+0.71269 

2 

— 

.17019 

- .56735 

- .61 300 

- .81113 

- .95736 

1 

+ 

7352 

+ .11108 

+ .15571 

+ .31108 

+ .57527 

6 

— 

572 

- 1160 

- 1935 

- 7563 

- .20752 

8 

+ 

27 

+ 71 

+ 116 

+ 1060 

+ 5010 

10 

— 

1 

3 

7 

102 

863 

12 





+ 7 

+ 111 

11 






11 

16 






+ 1 

Xi(l) 

+1.50038 

+1.62031 

+1.72298' 

+2.08616 

+2.51279 

Xi(0)* 

+ 

.37136 

+ .26571 

+ .19192 

+ .OII30 

+ .00121 
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TABLE 14 NUMERICAL VALUES OF THE COEFFICIENTS 
THE FUNCTION VALUES X^^ljr) AND X^ 


a^ r (7), AND 

(0;r) 


r 

- 7 2 = 3 

* 

5 

6 

-2 

+0.071289 

+0.095777 

+O.H9671 

+0.142474 

0 

+ .985721 

+ .974148 

+ .959391 

+ .941924 

2 

72766 

- 96430 

- .119654 

- .142393 

4 

+ 1840 

+ 3258 

+ 55067 

+ 7260 

6 

24 

57 

111 

190 

8 

+ 1 

+ 1 

+ 1 

+ 3 

XgU) 

+O.989062 

+0 .978117 

+O.964533 

+0 .949296 

*2(0) 

- .537431 

- .547702 

- .556050 

- .562257 


r 

-v 2 = 7 

8 

9 

10 

-2 

+0.163772 

+0.183325 

+0.201034 

+0.216893 

0 

+ .922384 

+ .901424 

+ .879642 

+ .857549 

2 

— .164669 

- .186536 

- .208085 

- .229432 

4 

+ 9837 

+ 12798 

+ 16148 

+ 19901 

6 

302 

- 450 

640 

880 

8 

+ 6 

+ 10 

+ 16 

+ 24 

X 2 (l) 

+0.933426 

+0.917893 

+0.903497 

+0.890893 

x 2 (o) 

- .566205 

- .567935 

- .567698 

- .565615 


r 

-7 

2 = 15 

20, 

25 

50 

100 

-2 

+0.27303 

+0.30282 

+0.31712 

+0.28828 

+O.22335 

0 

+ 

.75076 

+ .65357 

+ .55985 

+ .13410 

- .20363 

2 

— 

.33636 

- .44827 

- .56362 

- .91568 

- .79124 

4 

+ 

4532 

+ 8373 

+ .13671 

+ .50022 

+ .89250 

6 

— 

306 

768 

- 1596 

- .12202 

- .41949 

8 

+ 

13 

+ 43 

+ 112 

+ 1745 

+ .11556 

10 

— 

1 

2 

5 

166 

- 2129 

12 





+ 11 

+ 240 

14 





1 

21 

16 






+ 2 

x 2 (i) 

+0.86260 

+0 489086 

+0.96019 

+1.40297 

+1.81574 

x 2 ( 0 ) 

- 

.53563 

- .48567 

- .42431 

- .13884 

- .01312 
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TABLE 15 NUMERICAL VALUES OF THE COEFFICIENTS 
THE FUNCTION VALUES X^^ljr) AND dX^ 


a° (r), AND 

3,r 

(Ojr) /at 


X 

3 

4 

5 

6 


+0 .049768 

+0.065474 

+0 .080671 

+0.095329 


+ .996217 

+ .993384 

+ .989910 

+ .985786 


- 52769 

- 70252 

- 87661 

- .104980 


+ 1057 

+ 1877 

+ 2929 

+ 4212 


12 

27 

53 

92 

8 



+ 1 

+ 1 

I I 

+1 .000287 

+1.000066 

+0.999883 

+0.999742 

Eje|| 

-1.447434 

-1.427867 

+1.407450 

-1 . 386162 


"Nr 2 

r \ 

7 

8 

9 

10 

-2 

+0.109432 

+0.122969 

+0.135938 

+0.148343 

0 

+ .981091 

+ .975879 

+ .970201 

+ .964106 


- .122202 

- .139316 

- .156318 

- .173200 


+ 5725 

+ 7466 

+ 9433 

+ 11626 


145 

217 

- 308 

422 

8 

+ 2 

+ 4 

+ 7 

+ 10 

HM 

+0.999733 

+0.9999.13 

+I.OOO329 

+1.001021 

EES 

-1.364111 

-1.341380 

-1.318040 

-1.294172 



15 

20 

25 

50 

100 

-2 

+0.20229 

+0.24438 

+0.27699 

+0.35635 

+0.37464 

0 

+ .92878 

+ .88803 

+ .84416 

+ .61460 

+ .27353 


- .25562 

- .33423 

- .40827 

- .68495 

- .82807 


+ 2590 

+ 4548 

+ 6992 

+ .23918 

+ .56017 


141 

332 

641 

- 4395 

- .19704 

8 

+ 5 

+ 15 

+ 37 

+ 508 

+ 4388 

10 



2 

41 

681 

12 




+ 2 

+ 78 

14 





7 

BM 

+I.OO947 

+1.02683 

+1.05216 

+1.23184 

+1.53571 


-1.16949 

-1.04147 

- .91584 

- .42163 

- .07382 
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tabue 16.- udmosical taloes or ra cauricmrre r (r), aid the roicTio* vaiubs x^^ijr) 
AID xi^dir) aid a3$}ho i7 )/n 




+ 2349 

+ 67563 
+0 .994218 
- 68652 

+ 1097 

29 

+ 


+0 .999582 

+ .388021 


+ 9906 

+0.137045 
+ .975942 
- .136070 
+ 7544 

234 
+ 5 


+0.992656 
+ .403891 


+ 2240 

+0.19952 
+ .94636 
- .20137 
+ 1684 

79 

+ 2 


+O.98827 
+ .41853 


+ 4048 

+0.26279 
+ .90051 
- .26269 
+ 2956 

185 

+ 8 


+0.97338 
+ .43549 


+ 6318 

+0.32153 
+ .84069 
- .31828 
+ 4535 

357 

+ 18 
1 


+O.94971 
+ .45131 




+0.17284 

+ .45662 
+ .40843 
- .50209 
+ .16360 


6 

v + 2753 

8 

+ 291 

10 

21 

12 

+ 1 

14 


16 


V 15 

+0.82102 

14(0) 

+ .43598 


+0.19417 
+ .25147 

- .27937 

- .60518 
+ .599*1 

- .23640 

+ 5458 

849 

+ 96 

8 

+ 1 


+1.16944 
+ .14245 



10 

15 

+ 6543 

+ 1424 

+0.110949 

+0.16441 

+ .985522 

+ .96761 

- .112180 

- .16690 

+ 5292 

+ 1185 

143 

48 

+ 3 

+ 1 

+0 .998734 

+0.99668 

+1.788369 

+1.73792 



+ 2434 

+0.21549 
+ .94291 
- .22006 
+ 2093 

114 
+ 4 


+0.99393 

+ 1.68269 


+ 3632 

+0.26342 
+ .91191 
- .27117 
+ 32*5 

221 

+ 10 


+0 .99074 
+1.62227 


+0.10628 
+ .43712 
+ .6899* 

- .48639 
+ .12243 

- 1710 

+ 160 

16 


+0.98672 

+1.25100 


+0 .19798 

+ .49281 
+ .18816 
- .68181 
+ .39752 
- .11798 
+ 2220 
293 
+ 29 

2 


+1.11609 
+ .51419 
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TABLE IT.- NUMERICAL VALUES OF THE COKFFICIDnS a (7) AS WELL AS OF 

u,r 

THE I-FUNCTIOHS OF THE FIRST KHID FOR THE ARGUMENTS 1 AHD 0, 

WITH n = 6, 7, 8, 9, 10. 12 



+0.999935 
- .31737*+ 


10 

15 

20 

25 

+ 165 

+ 5001 

+ 94694 

+0.989966 
- 95268 

+ 3910 

91+ 

+ 3 

+ 57 

+ 1123 

+0 .14139 
+ .98042 
- .14252 
+• 879 

32 

+ 1 

+ 136 

+ 1979 

+0.18603 
,+ .96011 
- ,18796 
+ 1552 

75 

+ 2 

+ 268 

+ 3074 

+0 .22952 
+ .93789 
- .23256 
+ 2410 

145 

+ 6 

+0.999383 
- .3221+30 

+1.00132 
- .32872 

+0.99675 
- .33331 

+0 .99460 
- .33924 



n = 

6 

\ 2 



Nr7 

50 

100 

r 



-6 

+ 2977 

+0 .12226 

-4 

+0.11872 

+ .32925 

-2 

+ .40974 

+ .41829 

0 

+ .75313 

+ .12434 

2 

- .1+2313 

- .51561 

4 

+ 9124 

+ .28174 

6 

- 1122 

- 7755 

8 

+ 92 

+ 1354 

10 

5 

166 

12 


+ 15 

lb 


+ 1 

*6(1) 

+0.95890 

+0.80329 

Xg(0) 

- .3825*+ 

- .39599 



+ 879 

+ 5205 

+0.19542 
+ .50325 
+ .59 *+10 

- .52913 
+ .17076 

- 3192 
+ 1 + 01 + 

37 

+ 3 


+0.97926 

+1.99652 



+ 6383 

+0.25861 
+ .52093 
+ .38995 

- .58698 
+ .24863 

- .05770 

+ 880 

96 

+ 8 


n 

= 8 

25 

100 

+ 7 

+ 131 

+ 1816 

+0.17929 
+ .96481 
- .18035 
+ 1476 

72 

+ 2 

. 

+ 1719 
+ 7722 
+0.23556 
+ .51621 
+ .49454 

- .55330 
+ .20286 

- 4203 

+ 581 

59 

+ 5 

+0.99830 

+0.95851 

+ .28617 

+ .31+251 



+0.99965 
+ .2302U 
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ADDENDUM 


. NACA TM 1224 

/ 

LAME'S WAVE FUNCTIONS OF THE ELLIPSOID OF REVOLUTION 

By J. Meixner 

April 1949 


It has recently been brought to the attention of the NACA by 
Miss Gertrude Blanch of the Bureau of Standards, Department of Commerce 
that errors exist in the tabulated values appearing in tables 11 to IT 
of TM 1224. Miss Blanch notes that C. J. Bouwkamp, from whcm Meixner 
obtained the values presented, subsequently corrected them in tables- 
appearing in the Journal of Mathematics and Physics, vol. XXVI, no. 2, 
July 1947, pp. 88-91. 

In spite of the difference in symbols and notation in the two 
papers, reprints of tables I to IX included in the July 1947 issue of 
the Journal of Mathematics and Physics are attached for the use of those 
interested in receiving them. The NACA wishes to express its appreci- 
ation to the Journal of Mathematics and Physics for permitting these 
tables to be reproduced for this purpose. 
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TABLE I 


*« ■ 

Ao 

Ai 

At 

-10 

2.305040 

7.285254 

11.790394 

-9 

2.136732 

6.820888 

11.192939 

-8 

1.959207 

6.342739 

10.594773 

-7 

1.771184 

5.850492 

9.997253 

-6 

1.571156 

5.343904 

9.401958 

-5 

1.357357 

4.822809 

8.810735 

-4 

1.127734 

4.287129 

8.225713 

-3 

0.879934 

3.736870 

7.649318 

-2 ' 

0.611314 

3.172128 

7.084258 

-1 

0.319000 

2.593085 

6.533473 

0 

0 

2 

6 

1 

- 0.348602 

1.393206 

5.486800 

2 

- 0.729392 

0.773098 

4.996484 

3 

- 1.144328 

0.140119 

4.531027 

4 

- 1.594493 

- 0.505244 

4.091509 

5 

- 2.079934 

- 1.162478 

3.677958 

6 

- 2.599668 

- 1.831051 

3.289357 

7 

- 3.151841 

- 2.510421 

2.923796 

8 

- 3.733982 

- 3.200049 

2.578732 

9 

- 4.343293 

- 3.899400 

2.251269 

10 

- 4.976896 

- 4.607952 

1.938420 


TABLE II 


k * 

Ai 

Ai 

A* 

A« 

0 

12 

20 

30 

42 

1 

11.492121 

19.495277 

29.496855 

41.497757 

2 

10.990438 

18.994079 

28.995964 

40.997089 

3 

10.494513 

18.496395 

28.497321 

40.497988 

4 

10.003864 

18.002228 

28.000923 

40.000458 

5 

9.517981 

17.511597 

27.506763 

39.504497 

6 

9.036338 

17.024541 

27.014846 

39.010106 

7 

8.558395 

16.541109 

26.525161 

38.517282 

8 

8.083615 

16.061383 

26.037710 

38.026027 

9 

7.611465 

15 . 5 S 5448 

25.552488 

37.536339 

10 

7.141428 

15.113424 

25.069492 

37.048221 
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TABLE III 


QO 

Characteristic function X«(£) =■ ^ b*, Pin(t) 

o 



ba 

bt 

bt 

b . 

h 

bit 

-10 

0.944709 

- 0.728578 

0.110455 

- 0.007498 

0.00028 8 

- 0.000007 

-9 

0.951472 

- 0.684479 

0.094690 

- 0.005827 

0.000203 

- 0.000005 

-8 

0.958380 

- 0.635639 

0.079247 

- 0.004365 

0.000136 

- 0.000003 

-7 

0.965363 

- 0.581441 

0.064298 

- 0.003120 

0.000085 

- 0.000001 

-6 

0.972311 

- 0.521212 

0.050067 

- 0.002097 

0.000049 

- 0.000001 

-5 

0.979071 

- 0.454254 

0.036840 

- 0.001294 

0.000025 


-4 

0.985428 

- 0.379882 

0.024958 

- 0.000706 

0.000011 


-3 

0.991099 

- 0.297493 

0.014835 

- 0.000316 

0.000004 


-2 

0.995716 

- 0.206682 

0 . 006949 - 

- 0.000100 

0.000001 


-1 

0.998846 

- 0.107374 

0.001824 

- 0.000013 



0 

1.000000 






1 

0.998691 

0.114368 

0.001976 

0.000014 



2 

0.994509 

0.233927 

0.008138 

0.000118 

0.000001 


3 

0.987210 

0.356205 

0.018683 

0.000408 

0.000005 


4 

0.976790 

0.478301 

0.033563 

0.000979 

0.000016 


5 

0.963507 

0.597278 

0.052483 

0 . 001914 ' 

0.000039 


0 

0.947848 

0.710493 

0.074931 

0.003279 

0.000079 

0.000001 

7 

0.930440 

0.815971 

0.100266 

0.005114 

0.000143 

0.000003 

8 

0.911948 

0.912502 

0.127799 

0.007438 

0.000238 

0.000005 

9 

0.892980 

0.999612 

0.156881 

0.010250 

0.000369 

0.000008 

10 

0.874065 

1.077435 

0.186946 

0.013533 

0.000540 

0.000014 


TABLE IV 


QO 

Characteristic function &2nnJ > 2f l +i(£) 

o 


A * 

bi 

bt 

bt 

bi 

bt 

b ,, 

-10 

0.964429 

- 0.402104 

0.046184 

- 0.002528 

0.000081 

- 0.000002 

-9 

0.970923 

- 0.364430 

0.037096 

- 0.001858 

0.000054 

- 0.000001 

-8 

0.976877 

- 0.325710 

0.029954 

- 0.001312 

0.000034 

- 0.000001 

-7 

0.982232 

- 0.280082 

0.023016 

- 0.000882 

0.000020 


-6 

0.986936 

- 0.245730 

0.016934 

- 0.000556 

0.000011 


-5 

0.990948 

- 0.204851 

0.011752 

- 0.000322 

0.000005 


-4 

0.994236 

- 0.163056 

0.007499 

- 0.000104 

0.000002 


-3 

0.996784 

- 0.122359 

0.004197 

- 0.000069 

0.000001 


-2 

0.998586 

- 0.081179 

0.001852 

- 0.000020 



-1 

0.999051 

- 0.040320 

0.000159 

- 0.000002 



0 

1.000000 






1 

0.999604 

0.039016 

0.000447 

0.000002 



2 

0.998683 

0.078362 

0.001764 

0.000019 



3 

0.997105 

0.110098 

0.003902 

0.000063 

0.000001 


4 

0.994984 

0.152711 

0.006812 

0.000147 

0.000002 


5 

0.992373 

0.188112 

0.010436 

0.000281 

0.000004 


6 

0.989330 

0.222236 

0.014716 

0.000473 

0.000009 


7 

0.985910 

0.255039 

0.019595 

0.000733 

0 . 000016 - 


8 

0.982167 

0 . 2S0500 

0.025011 

0.001066 

0.000027 


9 

0.978150 

0.316012 

0.030908 

0.001477 

0.000042 

0.000001 

10 

0 . 97390S 

0.345386 

0.037230 

0.001909 

0 . 0000C2 

0.000001 
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TABLE V 


•0 

Characteristic function Xi(£) = buP*,(£) 

o 


k* 

A# 

b. 

tu 

b* 

b% 

An 

0 


1.000000 





i 

- 0.022875 

0.998525 

0.024445 

0.000206 

0.000001 


2 

- 0.046799 

0.993846 

0.048736 

0.000821 

0.000007 


3 

- 0.071286 

0.985722 

0.072766 

0.001840 

0.000024 

0.000001 

4 

- 0.095772 

0.974150 

0.096431 

0.003258 

0.000057 

0.000001 

5 

- 0.119671 

0.959391 

0.119654 

0.005067 

0.000110 

0.000002 

6 

- 0.142464 

0.941931 

0.142398 

0.007260 

0.000190 

0.000003 

7 

- 0.163759 

0.922394 

0.164677 

0.009837 

0.000302 

0.000006 

8 

- 0.183310 

0.901438 

0.186545 

0.012799 

0.000450 

0.000010 

9 

- 0.201017 

0.879661 

0.208098 

0.016150 

0.000640 

0.000016 

10 

- 0.216892 

0.857550 

0.229438 

0.019902 

0.000880 

0.000024 




TABLE 

VI 





Characteristic function A ' s ( S ) = ^ ta + i ^ fci + iU ) 

0 


k * 

b, 

4 * 

A* 

&7 

b% 

Ait 

0 


1.000000 





1 

- 0.016979 

0.999565 

0.017626 

0.000118 



2 

- 0.033587 

0.998287 

0.035224 

0.000470 

0.000003 


3 

- 0.049768 

0.996217 

0.052770 

0.001057 

0.000012 


4 

- 0.065475 

0.993406 

0.070263 

0.001877 

0.000027 


5 

- 0.080671 

0.989910 

0.087661 

0.002929 

0.000053 


6 

- 0.095328 

0.985786 

0.104979 

0.004212 

0.000092 

0.000001 

7 

- 0.109432 

0.981091 

0.122202 

0.005725 

0.000145 

0.000002 

8 

- 0.122970 

0.975879 

0.139316 

0.007465 

0.000217 

0.000004 

9 

- 0.135939 

0.970201 

0.156318 

0.009433 

0.000308 

0.000006 

10 

- 0.148343 

0.964106 

0.173200 

0.011620 

0.000122 

0.000010 


TABLE VII 


OO 

Characteristic function A' 4 (£) = ^ i) 

o 


*« 

Ai 

As 

b. 

b. 

A« 

An 

iu 

0 

1 

2 

0.000091 

0.000368 

- 0.013588 

- 0.027140 

1.000000 

0.999768 

0.999074 

0.013773 

0.027528 

0.000070 

0.000304 

0.000002 


3 

0.000834 

- 0.040653 

0.997918 

0.041266 

0.000084 

0.000006 


4 

0.001493 

- 0 . 05-1128 

0.990300 

0.054977 

0.001215 

0.000015 


5 

0.002348 

- 0.067563 

0.994218 

0.068051 

0.001897 

0 . 000029 


6 

0.003404 

- 0.080957 

0.991669 

0.082279 

0.002729 

0.000051 

0.000001 

7 

0.004663 

- 0.094312 

0.988048 

0.095855 

0.003712 

0 . 0000.80 

0.000001 

8 

0.006125 

- 0.107546 

0.984408 

0.109286 

0.001839 

0.000120 

0.000002 

9 

0.007806 

- 0.120900 

0.981162 

0.122811 

0.000122 

0.000171 

0.000003 

10 

0.009695 

- 0.134130 

0.976680 

0.130173 

0.007550 

0.000234 

0.000005 


■ ■ i 
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TABLE VIII 


oo 

Characteristic function X*ft) = b 2n +i Pzn+i ft) 

o 


k* 

i. 

b, 

b. 

b , 

b, 

ill 

bu 

0 

1 

2 

0.000068 

0.000273 

-0.011218 

-0.022423 

1.000000 

0.999854 

0.999418 

0.011294 

0.022584 

0.000046 

0.000213 

0.000001 


3 

0.000611 

-0.033609 

0.998690 

0.033863 

0.000478 

0.000004 


4 

0.001080 

-0.044772 

0.997673 

0.045127 

0.000850 

0.000009 


5 

0.001681 

-0.055906 

0.996367 

0.056373 

0.001327 

0.000018 


6 

0.002408 

-0.067006 

0.994772 

0.067596 

0.001910 

0.000031 


7 

0.003262 

-0.078066 

0.992888 

0.078793 

0.002598 

0.000050 

0.000001 

8 

0.004236 

-0.089080 

0.990719 

0.089958 

0.003391 

0.000073 

0.000001 

9 

0.005331 

-0.100043 

0.988263 

0.101088 

0.004290 

0.000105 

0.000002 

10 

0.006543 

-0.110949 

0.985522 

0.112180 

0.005292 

0.000143 

0.000003 


TABLE IX 


Characteristic function X«(£) = ^ P* n (£) 

o 


k 1 


bt 

b, 

bt 



bu 

bit 

0 

1 

2 

-0.000001 

0.000050 

0.000202 

-0.009535 

-0.019061 

1.000000 

0.999899 

0.999597 

0.009571 

0.019140 

0.000039 

0.000157 

0.000001 


3 

-0.000004 

0.000454 

-0.028580 

0.999094 

0.028702 

0.000353 

0.000003 


4 

-0.000010 

0.000805 

-0.038087 

0.998390 

0.038257 

0.000627 

0.000006 


5 

-0.000020 

0.001256 

-0.047578 

0.997486 

0.047802 

0.000980 

0.000012 


6 

-0.000035 

0.001807 

-0.057050 

0.996381 

0.057331 

0.001410 

0.000020 


7 

-0.000056 

0.002458 

-0.066501 

0.995076 

0.066845 

0.001919 

0.000032 


8 

-0.000084 

0.003207 

-0.075928 

0.993572 

0.076341 

0.002505 

0.000048 

0.000001 

9 

-0.000120 

0.004055 

-0.085326 

0.991867 

0.085816 

0.003169 

0.000068 

0.000001 

10 

—0.000165 

0.005001 

-0.094694 

0.989966 

0.095268 

0.003910 

0.000094 

0.000002 


- 900 
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